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This  appendix  is  intended  to  serve  both  as  a  source  of  immediately 
applicable  methodology  and  as  a  guide  to  the  underlying  gas  dynamic 
theory.  Those  interested  more  in  applying  the  methodology  than  in  der¬ 
ivations  and  comparisons  of  calculated  and  observed  results  will  find  the 
following  parts  of  this  appendix  of  particular  importance: 

Location  _ Description _ _ 

Section  IIA  The  simplest  and  fastest  method  of  estimating 

average  pressure  in  a  room  as  a  function  of 
time;  adequate  for  many  purposes. 

Section  IIB3  Two  different  methods  for  a  step-by-step,  hand 
calculation  providing  average  room  pressure  as 
well  as  dynamic  pressure  in  the  opening,  valid 
for  all  flows  through  a  single  opening  into  a 
single  room  when  the  outside  pressure  is  known 
as  a  function  of  time.  Either  method  may  be 
used  but  Method  F  is  simpler  for  inflow. 

Method  D  for  outflow;  they  may  be  combined  in 
use. 

Section  IIC  Formulas  for  geometrical  extent  of  the  jet 

created  inside  the  room  by  inflowing  air  and 
for  dynamic  pressure  distribution  within  it. 

Table  E-3  A  computer  program  to  calculate  average  pres¬ 

sure  within  a  single  room  and  dynamic  pressure 
in  as  many  as  eight  openings  as  functions  of 
time  when  a  single  room  has  openings  into  sev¬ 
eral  different  pressure  fields  (e.g.,  a  room 
with  front,  rear,  and  side  windows  struck  by 
a  blast  on  the  front  wall). 

Section  IVB  A  step-by-step  hand  calculation  of  average 

room  pressure  and  dynamic  pressure  in  each 
opening  during  filling  of  a  single  room 
through  one  or  two  openings  into  separate 
pressure  fields. 
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When  using  the  methodologies  noted  above,  the  user  may  find  the  mean¬ 
ings  of  symbols  under  Notation  at  the  end  of  this  appendix.  Subscripts 
not  explained  there  refer  to  physical  spaces,  i.e.,  subscript  1  indicates 
quantity  is  measured  outside  the  room;  other  odd-numbered  subscripts  refer 
to  interior  of  rooms  and  even-numbered  subscripts  refer  to  connecting  ducts 
or  openings. 
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Appendix  E 


I  Introduction 


ROOM  FILLING  FROM  AIR  BLAST 
By  J.  R.  Rempel 


Mien  a  blast  wave  strikes  a  building,  even  should  the  structure  with¬ 
stand  the  initial  impact,  the  resulting  inflow  of  air  through  windows  and 
other  openings  can  be  critical  in  determining  the  safety  of  any  people 
sheltered  by  the  structure  and  in  determining  the  response  of  the  struc¬ 
ture  itself  to  the  blast  impact.  Although  the  physical  laws  obeyed  by 
moving  gases  are  well  known  and  the  course  of  the  inflow  in  filling  the 
building  can  in  principle  be  calculated  completely,  any  such  calculation 
is  far  too  lengthy  to  be  practical  for  most  purposes;  fortunately,  sim¬ 
plifications  can  be  introduced  which  greatly  shorten  the  labor  of  estimat¬ 
ing  effects  of  the  blast  inside  the  building  and  which  give  results  in 
good  to  fair  agreement  with  experiments  done  with  small  models.  In  gen¬ 
eral,  the  effect  of  the  inflow  is  to  provide  a  stream  of  fast  moving  air 
in  the  shelter  space  which  may  (1)  endanger  shelterees  by  hurling  them 
against  large  relatively  fixed  objects  or  by  hurling  objects  against  them,  . 
and  (2)  provide  a  back  pressure  on  the  inner  surfaces  of  structure  walls 
countering  the  blast  pressures  on  their  outer  surfaces. 

Several  factors  enter  into  the  calculation:  the  pressure  outside 
each  wall  with  openings  and  the  time  each  opening  becomes  available,  the 
area  occupied  by  each  opening  and  the  volume  of  each  room,  the  number  of 
connected  rooms  and  the  area  of  each  connection,  and  the  ambient  pressure 
and  temperature  in  the  building  before  the  blast  strikes.  Perhaps  the  first 
of  these  to  consider  is  what  proportion  of  the  wall  exposed  to  the  blast 
is  open.  If  this  fraction  is  greater  than  one  half,  the  shock  front  lead¬ 
ing  the  blast  wave  will  pass  into  the  building  only  slightly  weakened  and 
subsequent  inside  pressure  should  be  estimated  from  a  knowledge  of  shock 
position  and  the  laws  of  shock  reflection.  Methods  appropriate  to  this 
case  are  only  touched  upon  here.  On  the  other  hand,  should  the  fraction 
be  less  than  one  tenth,  clearly  the  filling  is  not  a  shock  process  and 
the  methods  treated  here  are  quite  pertinent.  Unfortunately,  in  many 
applications  the  fraction  of  open  area  will  lie  between  these  two  ex¬ 
tremes  and,  in  these  cases  after  the  room  filling  calculation  has  been 
completed  according  to  the  methods  suggested  here,  some  thought  must  be 
given  independently  to  the  influence  of  the  entering  shock  fi'ont. 
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Wlion  the  source  of  the  blast  wave  is  an  explosion  and  the  location  of 
the  building  in  relation  to  the  point  of  explosion  is  known  or  postulated, 
"free  field"  pressure  histories  at  the  building  site  can  bo  found  in  stan¬ 
dard  references  ,^'^'^^and  from  those  histories  well  known  methods^ are 
available  to  derive  approximate  histories  on  the  outside  of  the  walls  of 
the  building.  Briefly,  these  methods  account  for  a  short-lived  peak  of 
pi'essure  created. by  the  impact  of  the  front  upon  the  wall  nearest  the  ex¬ 
plosion,  the  relatively  fast  erosion  of  this  high  pressure  to  a  level  which 
is  the  sum  of  the  free  field  pressure  plus  a  drag  pressure  on  the  wall 
due  to  the  high  winds  behind  the  blast  front ,  This  quasi— steady  pressure 
then  decays  slowly  to  zero  as  the  blast  wave  moves  onward  past  the  struc¬ 
ture  , 

Ordinary  window  glass  breaks  rather  quickly,  i  ,e ,  ,  within  8  ms  (milli* 
seconds)  or  less  when  struck  by  blast  overpressure  of  1  psi  (pound  per 
square  inch)  or  more  Doors  may  withstand  outside  pressure  longer,  or 
even  altogether.  The  time  an  opening  becomes  available  with  respect  to 
the  first  impact  of  the  blast  upon  the  building  becomes,  then,  the  break¬ 
ing  time  plus  the  time  required  by  the  wave  to  travel  from  the  wall  near¬ 
est  the  explosion  to  the  opening.  If  the  opening  is  in  the  wall  nearest 
the  explosion,  travel  time  is  of  course  zero.  Strictly,  the  decay  of  the 
blast  wave  overpressure  which  occurs  during  this  time  must  bo  taken  into 
account,  but  when  the  blast  arises  from  a  nuclear  explosion  of  yield 
greater  than  a  few  kt  (kilotons)  this  decay  is  slight  and  negligible;  that 
is,  a  single  "free  field"  pressure  history  for  all  openings  may  be  assumed. 

It  should  be  emphasized  that  the  methods  given  here  are  simplified 
and  their  use  can  lead  only  to  estimates.  They  are  intended  to  provide: 

(1)  calculations  applicable  to  hand  computation  by  those  untrained  in  gas 
dynamics  and  (2)  approximate  results  useful  until  more  careful  calcula¬ 
tions  are  made.  Only  in  the  case  of  the  simplest  structural  configura¬ 
tions  and  the  simplest  pressure  history  shapes  can  limits  of  error  be  sug¬ 
gested  for  these  results.  Such  cases  will  be  the  subject  of  the  discus¬ 
sion  immediately  below. 


II  Classical  Nuclear  Blast  Wave  Incident  upon  a  Single  Room 
A.  Estimation  of  Inside  Pressure  History 

The  "classical"  blast  wave  from  nuclear  explosions  consists  of  a  steep 
pressure  front  or  rise  followed  by  a  long-lasting  decay  phase  during  which 
the  pressure  in  the  wave  falls  to  zero.  It  is  accompanied  by  high  winds 
giving  rise  to  dynamic  pressure  against  objects  in  the  stream.  Striking 
a  wall  at  normal  or  near  normal  incidence,  it  creates  a  high  pressure 


References  are  listed  at  the  end  of  this  Appendix. 
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zone  at  the  surface,  which  however  is  rapidly  eroded  as  relief  waves  move 
across  the  wall  face  from  the  edges.  As  a  first  approximation  in  the  cal¬ 
culation  of  room  filling,  this  reflected  phase  can  usually  be  neglected 
Following  the  decay  of  the  high  reflected  pressure,  the  quasi-steady  pres¬ 
sure  (free  field  plus  dynamic)  remains  against  the  wall  for  hundreds  of 
milliseconds  to  several  seconds,  depending  on  explosive  yield.  Generally 
filling  is  complete  before  this  quasi-steady  pressure  has  fallen  more  than 
a  few  percent;  hence,  as  a  first  approximation  the  outside  pressure  may 
often  be  considered  constant  and,  if  only  a  single  wall  with  opening  is 
exposed  to  the  blast,  the  time  AT  (in  milliseconds)  to  complete  filling 
may  be  computed  as  the  ratio  where  V  is  room  volume  in  cubic  feet  and 
A  IS  area  of  opening  in  square  feet.t  The  average  room  pressure  at  any 
time  t  during  the  filling  process  is  then  simply  the  fraction  of  the  quasi¬ 
steady  outside  pressure  given  by  the  ratio  For  the  purposes  of  this 

calculation,  areas  of  several  openings  in  the  same  wall  should  be  added 
together  to  form  the  quantity  A** 


In  case  there  are  two  or  more  walls  with  openings  exposed  to  the  blast 
and  each  such  wall  sustains  a  different  outside  pressure  history  (as  will 
happen,  for  example,  when  the  drag  coefficient  is  different  for  two  walls) 
the  calculation  is  more  complicated  but  first  estimates  of  filling  time 
and  average  inside  pressure  during  filling  may  be  found  by  adding  interior 
pressures  calculated  as  if  each  wall  alone  were  exposed.  As  an  example 
consider  a  room  of  volume  30’xl0*xl0'  =  3000  ft^  in  which  the  front  wall 
has  total  openings  of  36  ft^  and  side  walls  have  total  openings  of 
60  fts.  The  ratios  for  the  front  and  side  walls  are  41.7  and  25  ft, 
respectively.  If  the  quasi-steady  overpressure  on  the  front  wall  is  lo'psig 
(pounds  per  square  inch  -  gauge)  and  on  the  side  wall  8  psig  and  if  further 
the  side  wall  opening  becomes  available  10  ms  after  the  first  blast  im¬ 
pact,  then  the  average  inside  pressure  will  be  approximately  as  shown  in 
Figure  E-1  by  the  heavy  line  OAFG.  In  other  words  the  room  will  fill  in 
approximately  24  ms.  Lines  OAC  and  DE  represent  filling  rates  through 
front  and  side  walls, respectively;  and  ordinates  of  OAC  and  DE  are  added 
to  form  the  line  OAFG.  Of  course  after  the  average  inside  pressure  ex¬ 
ceeds  8  psi  there  will  be  outflow  through  the  side  wall;  to  allow  for 
this  loss  the  line  FG  has  been  placed  between  the  outside  pressure  at 
the  side  wall  (8  psig)  and  the  outside  pressure  at  the  front  wall  (10  psig). 
The  ordinate  at  FG  is  closer  to  8  psig  than  to  10  psig  because  the  area 
of  the  opening  in  the  side  wall  is  greater  than  that  in  the  front  wall. 

The  lino  FG  is  intended  to  represent  the  final  quasi-oqulllbrlum  pro.ssuro 
In  the  room. 


*  The  experimental  justification  of  most  of  the  procedures  described  in 
this  section  is  demonstrated  later  in  Figures  E-4,  E-5,  E-6  and  E-8. 
t  Empirical  relationship,  dimensionally  inconsistent.  Meanings  of 
symbols  as  used  in  this  Appendix  are  defined  as  introduced  and  under 
Notation  at  end  of  Appendix. 
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_  _  OUTSIDE  PRESSURE  (FRONT) 

X  ^  final  equilibrium 

®  PRESSURE  IN  ROOM 

_ OUTSIDE  PRESSURE  (SIDE) 


TIME — ms 

P.3UR.  R-1  APPBOX,^ATE  P.LLING  BATES  THROUGH  TWO  WAULS 

These  simple  calouletions  "°^"°'XwhL*tte  wave  is  non-elasslcal 
asts  present.  Under  these  eondltiens  the 

.ore  retailed  methods  set  forth  below  must  be  followed. 


B  ^  General  Case 

justification  for  use  of  the  simple  Ttfresrits 

xperience  with  a  step-by-step  prLciples  of  steady  isen- 

ith  experiments.  This  calcula  ion  p  intei-vals.  Conditions  com- 

ropic  flow  in  ducts  in  access  ve  small  the  next 

)uted  for  the  end  of  one  time  s  ep  along  with  the  assumption 

,tep.  Conservation  of  ene  gy,  constant  specific  heats  deter- 

;hat  the  air  behaves  as  a  per  temperature  and  density,  as 

nine  the  thermodynamic  vaiiablo  ,  l  >  expressions  lending 

veil  as  the  wind  speed  through  for  the  laws  of  conserva- 

themselves  to  simple  calculation  c  leading  to  somewhat  dif- 

ripp  Of  ener.v  Atr^r^rr^rrs  re?T  upoh  ceruiin 

ferent  results  1,.,  3nd  these  approximations  usually 

approximations  to  the  eomparison  with  whioh  the  approxima- 

introduce  errors  into  the  res  , ,^„troov  perfect  gas  behavior  and 

tions  arising  in  the  assumptions  of  isentropy.  pert 

constant  specific  heats  are  negligible. 
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1.  Inflow.  Figure  E-2  shows  the  icicnlizccl  room'  with  a  single  open¬ 
ing  struck  head-on  by  a  blast  wave.  Three  regions  are  noted:  the  outside 
(p,  tl^  doorway  ©which  serves  as  a  duct  connecting  the  outside  with  the 
room©.  In  order  to  make  the  calculations  tractable,  uniformity  of  con- 
ditions  in  each  of  the  two  regions (l)  and  ©  and  over  the  cross  section  of 
region©  IS  assumed;  furthermore,  during  each  small  time  interval  At,  steady 
conditions  arc  assumed  in  each  region.  During  the  aforementioned  quasi¬ 
steady  state  outside  the  building  these  assumptions  are  probably  valid  for 
region  ©but  they  clearly  introduce  error  if  the  reflection  or  diffraction 
phase  lasts  an  appreciable  time,  for  during  that  episode  relief  waves  are 
moving  into  the  region  from  the  edges  of  the  building  as  well  as  from  the 
doorway  itself  causing  rapid  fluctuations  in  wind  speed  and  pressure. 

(Some  account  is  taken  of  changes  in  pressure  during  the  diffraction  epi¬ 
sode  by  the  standard  techniques  of  estimating  outside  pressure.)  Similar 
remarks  c.an  be  made  concerning  regions  @  and  ©,  but  if  we  are  content  to 
deal  with  average"  pressure  and  speed  in  those  two  regions,  we  may  apply 
the  step-by-step  isentropic  analysis,  However,  our  present  methods  do 
not  provide  for  any  apportionment  of  gaseous  energy  in  region  ©  between 
streaming  kinetic  energy  and  internal  energy;  for  simplicity  of  calcula¬ 
tion  it  will  be  treated  as  entirely  internal  at  all  times,  which  will 
cause  overestimation  of  pressure  and  neglect  of  winds  within  the  chamber 
In  evaluating  the  wind  threat,  the  speed  and  dynamic  pressure  in  the 
duct  ©  must  be  regarded  as  the  upper  bounds  on  wind  speed  and  dynamic 
pressure  in  region©.  Later,  methods  will  be  given  for  estimating  change 
in  dynamic  pressure  as  the  wind  moves  into  the  room. 


FIGURE  E-2  THE  FILLING  CHAMBER 
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In  writing  the  conservation  equations,  two  views  can  be  taken  of  con¬ 
ditions  in  region®.  On  the  one  hand  pressure,  density  and 
be  those  of  the  free  field  behind  the  blast  front  or,  on  the  other  han J , 
the  air  upstream  of  the  opening  may  be  treated  as  stagnate  at  a  Pressure 
above  free  field  either  by  the  amount  of  the  reflected  pressure  or  by  the 
amount  of  the  product  of  the  drag  coefficient  and  dynamic 

vided  drag  coefficients  are  known,  the  second  view  is  more  simply  applied, 
especially  when  the  blast  front  does  not  meet  the  wall  head-on .In  what 
follows, the  pressure,  ,  and  density,  pi  ,  in  region®  will  be  those  o 
stagnate  air  outside  the  wall.  The  work  done  in  moving  a  mass  element 
in  region®  through  a  small  distance  Ax  toward  region®  is; 


P  A  dx  =  P  AV  =  —  •A"' 
11  1  1  Pi 


„hor«  A  is  the  cross  sectional  area,  and  AV  Is  the  volume  occupied  by 
Am  in  resion®.  The  mass  element  carries  with  it  the  internal  energy  i 


had  in  region®,  i.e., 


1  .  .  Am, 

Y  -  1 

Where  Y  la  the  ratio  ol  specllic  heat  at  constant  pressure  to  the  specillc 
heat  a?  constant  volume.  (The  perfect  gas  equation  ol  state  is  assumed; 
see  Ref  5  )  If  the  flow  into  the  room  is  steady,  energy  conservatio 
requires  that  the  same  total  energy,  specifically,  the  sum 


P  P  *1 

1  .  1  .  Am  .  Am  =  Y .  _i 

y  -  1  T  p^  Y--  1  Pi 


^  J.  .  Am 


be  given  up  wJ-an«gion  (^  during 

::rarrrrgror®  equal  The  worh  done  in  region©  in  pressing  the  mass 

element  toward  region®  is 


'  2  Am 


the  internal  energy  in  the  element  is 
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Y  -  1  P2 

whei'e  the  subsci'ipt  Q)  denotes  conditions  in  region  (^.  Since  the  air  is 
flowing  into  the  room,  however,  the  element  in(2)also  carries  (streaming) 
kinetic  energy  of  amount 

1  2 

where  u  designates  particle  or  material  speed.  Thus,  if  conditions  are 
not  changing  too  fast,  we  can  write  (cancelling  out  the  factor  Am): 


2  12 
—  +  —  ^ 
Po  2  2 


To  apply  conservation  of  linear' momentum  we  consider  a  control  sur¬ 
face,  shown  dashed  in  Figure  E-3 ,  Neglecting  gravitational  forces,  the 
x-component  of  the  force  integrated  over  this  surface  must  equal  the  rate 
of  flow  of  the  x-component  of  momentum  out  of  the  volume  plus  the  rate  of 
increase  of  x-momentum  within  the  volume.®  Neglecting  frictional  and 
viscous  effects,  the  only  force  on  the  surface  is  the  thermodynamic  pres- 


FIGURE  E-3  FIRST  CONTROL  SURFACE 
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sure  exerted  normal  to  the  surface  in  an  inward  direction.  Across  the 

ZoU  of  the  duct  this  pressure  is  taken  to  be  everywhere  else  on 

the  surface  it  is  approximately  P^ .  (Around  the  entrance  to 

.re.riure  on  the  wal]  wi  ]  ]  be  less  than  Pi,  but  the  mass  flow 

tst  sot  to  Ta-  hi;<hly  sensitive  to  corrections  made  for  thi  s  ef  fee  J , 

in  order  to  simplify  calculation,  uniform  pressure  Pi  will  oe 

The  mass  speed  uo  however  would  be' reduced  by  this  correction.)  s. 

r.  auasi^rt^dy  episode  .o«entu«  within  the  surface  will  chan.o  only 

‘“w^y  aid  we  aLuie  that  it  is  in  fact  constant  during  each  ™ 

tl.e  Lcrement  At.  (This  assumption  is  clearly  false  ‘he  shoch 

diffraction  episode  due  to  the  presence  of  fast  moving  reflected  shock 

iioiii Tirnih  the  surface.  During  the  guasi-steady  episode  the  assump- 

tirierficts  the  relief  wave  spreading  into  the  stagnate  air  -tside  the 

opening.)  Thus,  carrying  out  the  integration  over  the  sur  ace, 


P  COS 


6  da  =  (P.  -  P,) 


T  Q  •  e  flTip-le  betwcen  the  inward  normal  to  the  surface 

in  this  integral  6  is  on  the  control 

and  the  positive  k-direction  and  da  ^  through 

surface.  The  total  momentum  in  the  positive  x 

the  surface  is  simply 

pg  A  "'g 


Hence,  momentum  conservation  reduces  to 


/ 

p  _  p  =  u„ 

2  ^2  2 

Plnally  we  note  that  even  in  f-  PTesenee  of  "---i/n-^rrairlough 
Shocks  the  isentropic  equation  of  state  of  a  periec  g 
for  this  approximate  calculation;  hence 


Po  =  P 


P  1  1/Y 


1  P 


(See  Ref  .  7)  . 

VI  /IN  /vN  and  f3)  mav  be  solved  for  P2  >  P2 ' 

Given  Pi  and  Pi ,  Eqs.  (1),  (2),  and  may 

and  U2.*  The  result  can  be  written  as; 


.K  For  diatomic  gases  like  air,  y=lA\  see  Ref.  8. 
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I 


Y  +  1 


Y  +  1 


2  Y  -  1 

which  is  independent  of  p  .  .  If  y  =  — A  =  - - 

P  Y  +  ^ 


and  B  = 


_  — £X_  Eq.  (4)  can  be  put  in  the  form 
Y  +  1 


By  ^  =  y  +  A  (5) 

When  A  and  B  have  the  values  stated  above,  Eq.  (5)  has  two  solutions, 
one  of  which  is  y  =  1  and  the  other  is  y  =  0.1912.  The  second  solution 
is  the  only  one  of  interest  here  and  will  be  designated  y  . 


To  continue  the  calculation  must  be  known.  This  value  can  be 
found  from  the  Rankine-Hugoniot  relations  and  knowledge  of  the  strength 
and  angle  of  incidence  of  the  original  shock  front  (Ref.  1,2),  or  it  can, 
with  enough  accuracy  for  incident  shock  strengths  less  than  15  psi,  be 
computed  from  ambient  conditions  using  the  isentropic  equation  of  state , 


1  .e . , 


Pi  =  Po 


P  ll/Y 


where  P  and  p  are  ambient  pressure  and  density,  respectively.  With  p^ 
known,  air  density  and  pressure  in  the  opening  can  be  calculated  from; 

p^=p^y^l/V  (7> 


=  y  P-. 

2  o  1 


from  which  wind  speed  becomes: 


P,  -  2 

1  2  ^ 


-1  -  >'»> 


The  mass  flow  into  the  room  (3)  can  be  written 
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^  =  P2  ''2  ^2^^ 


=  P  (1  -  y  )  y  p 
1  o  o  "^1 


A2  At 


If  V3  is  the  volume  of  the  room  and  the  prime  is  used  to  denote  conditions 
in  the  room  at  the  beginning  of  time  step  At,  then  average  density  in  the 
room  p„  at  the  end  of  At  can  be  written  as: 

O 


/  Am 
•^3  "  ^3  V„ 


To  find  pressure  in  the  room  at  the  end  of  the  time  step  we  assume  that 
all  the  energy  lost  in  region  appears  as  an  increase  of  internal  en- 
ergy  of  the  gas  in  the  room,  i.e., 


— I —  .  —  .  ^  — 

Y  -  1  Y 


P  >  p 

3  3 


which  can  be  solved  to  give  in  terms  of  known  quantities: 


Am  / 

P  =  -  •  —  +  P^ 

^  Pi  ^3  " 


(12) 


At  any  time  air  temperature  within  the  room  can  be  calculated  from 
the  perfect  gas  law: 


^3  = 


is  the  gas  constant^  for  air  in  the  appropriate  units  (e.g.,  in 
metric  units  R  =  0.3028  joule/g^C) .  The  quantity  T3  can  reach  high  values 
as  a  result  of  the  compression  existing  behind  the  shock  and  within  the 
roomj  however,  if  airflow  to  the  outside  (outflow)  is  maintained,  the 
relaxation  of  pressure  following  the  passage  of  the  front  will  return 
room  gas  temperatures  to  safe  levels  before  injury  to  occupants  is  likely. 
Only  the  long  lasting  increase  in  room  temperature  resulting  from  fires 
will  normally  be  a  threat  to  the  shelterees. 
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Conditions  existing  in  the  chamber  at  thr>  •  • 

;.ro  used  in  tlio  lorosolng  calculations  only  i„ 

not  inllucnco  duct  parameters  because  tranLcnti  ■  ' 

consideration.  Transient  phenomena,  for  er,m,"le  d''t 

Of  flow;  that  is,  if  P  >  p  .  oxamplo,  determine  the  direction 

o...er.lse  Horn  Is  outwa^:.  \ “irn:::::: 

the  room  leads  to  an  accumulation  of  error  in  *  “‘•"Wnatlne  from 

pressure  as  can  be  soon  from  comparisons  between  00^0^11“” 

rr:\i"i:rr“ 

Chamber  is  shown  in  an  inset^in  the' tl^rr"^;rtL“r 

tion  produces  Curve  F  of  tho  fimi,  "  foregoing  calcula- 

e-6  labelled  ■'owtern.al  Lst.'T;  '  ; 

by  means  of  a  pitot  tube  oriented  wl  urrespect'’'trjhe*"t  "'r' 

\vi til  tho  orientation  of  -Hirk  j  to  confonn 

in  Figure  E-4  is  side-on  overprLsLo'\r'L"'°"''\''°°'”’ 

In  each  case  the  calculation  IrUf  unobstructed  shock  tube,^ 

observation  but  eventually  sh  ^  yields  pressure  in  agreement  with 

Of  measurement!  .uL“h  !he  r"- o-ess 
tude  of  At  for  tL“  eaL!!»t  ™  O'-  1---  M.gnl- 

sound  signal  across  the  room.'””*  *”*  one-quarter  the  transit  time  of  a 
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OVERPRESSURE  (PSI) 


FIGURE  E-4  FILL  HISTORIES  FOR  SIDE-ON  INCIDENCE 


'V 

u 
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OVERPRESSURE  (PSI  ) 


figure  e-5  fill  histories  for  a  face-on  model 
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FIGURE  E-6  FILL  HISTORIES  FOR  A  REAR  FILL  MODEL 
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Choice  of  the  size  of  At  is  somewhat  arbitrary  except  that  (1)  values 
much  greater  than  sound  transit  time  will  give  a  false  idea  of  the  degree 
of  irregularity  in  the  fill  process  and  (2)  enough  steps  should  be  taken 
to  make  it  possible  for  the  influence  of  variations  in  to  be  shown  in 
the  results.  The  length  of  the  bar  labelled  in  Figures  E~4 ,  E-5,  and 
E-6  represents  the  sound  transit  time  across  the  longest  room  dimension. 

The  existence  of  significant  theoretical  errors  in  our  treatment  of 
flow  into  a  room  by  quasi-steady  analysis  is  clearly  revealed  by  consider¬ 
ing  the  single  control  surface  formed  by  superposition  of  that  shown  in 
Figure  E-3  and  that  indicated  with  dashed  lines  in  Figure  E-7.  Such  a 
surface  coincides  with. the  inner  surfaces  of  the  room  and  passage  and 
extends  into  quiescent  air  outside.  Under  our  hypotheses  there  is  no 
flow  through  this  surface  anywhere  and  no  change  of  momentum  within  it, 
yet  the  surface  integral  of  the  x-component  of  pressure  over  the  bound¬ 
aries  does  not  vanish. 


FIGURE  E-7  SECOND  CONTROL  SURFACE 
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This  absurdity  can  be  avoided  in  one  or  both  of  two  ways.  A  term  may 
be  added  to  the  right  side  of  Eq.  (2)  to  account  for  the  changing  flow 
pattern  within  the  control  surface  shown  in  Figure  E-3  or  a  term  may  bo 
added  to  the  left  side  of  Eq.  (2)  to  account  for  the  possible  nonuniform- 
^ -j-y  pi*essure  over  the  boundaries  of  the  surface.  As  the  room  begins 

a  rarefaction  wave  moves  back  into  the  high  pressure  gas  outside 
the  first  doorway  bringing  more  and  more  gas  into  motion  toward  the  open¬ 
ing.  In  other  words,  the  neglect  of  the  rate  of  change  with  time  of  the 
momentum  within  the  control  surface  outside  the  room  may  bo  at  least  one 
cause  of  the  contradiction  noted  above.  Any  attempt  to  calculate  a  cor¬ 
rection  for  this  effect  would  certainly  add  to  the  complexity  of  those 
simplified  procedures;  furthermore,  the  degree  of  agreement  between  ob¬ 
servations  and  theory  of  Method  F  (shown  in  Figures  E-4,  E-5,  and  E-6) 
suggests  that  the  added  effort  to  account  for  the  rarefaction  wave  may 
not  be  needed  to  achieve  the  desired  degree  of  accuracy. 

Some  writers^®' use  the  equation 

P  =  P'  (13) 

2  3 

instead  of  Eq.  (2).  Justification  for  Eq.  (13)  is  based  on  analogy  ^1'^^ 
the  treatment  of  flow  into  a  large  chamber  steadily  being  evacuated. 
Equation  (13)  of  course  provides  continual  coupling  between  flow'  condi¬ 
tions  and  conditions  in  the  room. 


Using  Eq.  (13)  we  derive  the  pressure  buildup  inside  the  room  in 

the  following  way.  Substituting  Eq.  (3)  into  Eq.  (1),  then  replacing 

P„  with  Po  according  to  Eq.  (13)  and  solving  the  resulting  equation  for 
2 


U2,  we  find: 


^^2 


Y-1  P 


1/2 


1  \  \  1 

from  which  we  calculate  the  mass  inflow  in  the  time  increment  At; 


=  K  p2  u^  At 


=  K  p 


2y 


1\Y-1  P. 


1-  U] 


*2  “ 


(14) 


Whenever  Eq.  (13)  is  employed,  empirical  corrections  must  be  made  to  re¬ 
duce  the  calculated  inflow  rate  or  the  calculation  will  not  yield  realis' 
tic  values  for  room-pressure  in  small  experimental  models;  the  simplest 
correction  is  the  discharge  coefficient,^^  represented  by  the  factor  K 
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in  Eq  (14).  Investigators  at  the  IIT  Ko.-;oaivh  lu.-a  i  ux/' '  i.-uua 

it  necessary  to  use  the  value  K=0.7  to  reconcilo  corapuUHl  pros.sui-o  i  isoh 
in  small  models  with  those  measured;  Curves  D  in  Mijures 

bLnLuced  by  a  calculation  based  on  Eo.  (13),  with  K.O.l.durins 
inflow  and  K  =  1.00  during  outflow. 

The  value  of  the  discharge  coefficient  is  usually  discussed  in  con¬ 
nection  with  boundary  layer  thickness  and  the  Reynolds  number. 
sJoild  bfnoted  therefore  that  the  relatively  good  agreement  between 
Se  oLerved  room  pressures  and  Curves  D  was  obtained  in  very  small 
loLls  and  not  In  lull-sited  rocs,  lo  provide  an  estimate  of  the  in¬ 
fluence  of  the  value  of  K  on  calculated  pressuro  rise,  Curve  ,  aso 
fn  t^value  K  =  1.00  durlne  both  outflo.  and  inflow,  has  been  entered 
"  rraurf  E!!!  the  flow  Into,  lull-sited  rooms,  presumably,  the  Eey- 

nolds^umber  will  be  larger  and  the  discharge  coefficient  more  nearly 
equal  to  1.00  than  in  the  flow  into  small  models. 

Finally,  the  pressure  increment  during  the  interval  At  is  found  by 
substitution  of  Eq.  (14)  into  Eq.  (12): 


^3  "  ^3  ~  ^  ^  ^1  V-l  ’  Vl/  J 

Mellchar'^.'*  omits  both  the  factors  K  and  y  before  the  right-hand 
=  Of  EC  (15)  Which  is  equivalent  numerically  to  making  K-0.7  and 
V  1  4  He  attempts  to  justify  this  procedure  on  theoretical  grounds 
L-ccltS  :ith  Lndary  layer  theory  -  Some 

“aLTi;  Trl^/kL  ;f  sound  across  the  room. 

Pursuing  the  analogy  between  room  filling  and  steady  flow  into  a 

wMch^tsenLopL''flow\nto  tie  room  achieves  the  maximum  mass 
for  which  the  flow  speed  equals  local  sound  speed.  Therefore,  to 
risTritical  ratio,  we  can  differentiate  Eq  (14  with  -pe-  ^ 

set  the  result  equal  to  zero  and  solve  for  (Ps^^l^crit. ’ 


p;\^‘y  12  ^2 


1  /  crit. 
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TLTJZT.nl  S  t"t:inorriursU.uUni%.c  cnuinl  r„tlo 

into  Eq.  (14),  i.e.:  y+i  l 


r 


(AnO. 


choked 


=  K 


Y  Pi  Pi 


2 


A  At 


Because  the  mass  IbeylnrEq.  <W)^ls  called 

CcLred^e4:h::r 

equals  0.5283;  thus,  assuming  ambi  P 
outside  pressure  is 


_  27.83  psia  or  13.13  psig 

C  oo  o 


0.5283 

peak  overpressure  In  °i  '“’’""“orJccordlnr^  the'^oregoing: 

T?-'^  and  E-6  rose  above  13. lo  psig, 

the;ry,  choked  How  should  not  have  occurred. 

A  similar  degree  rtimlnritet  “  ^7  It^  model 

is  found  in  the  results  s  ow  except  that  the  transient  fluctu- 


p  =  13.58  psia 

°  ,  ‘ 
0  =  0.0672  Ib/ft 

^o 

Y  =  1.4 

X  =  0.821  ft^ 

3 

V  =  27.0  ft 
3 


ambient  pressure 
ambient  density 
ratio  of  specific  heats 
area  of  opening 
volume  of  model  room 


Phe  traction  ot  the  Impacted  wall  area  occupied  hy  the  opening  Is  slightly 
less  than  one-tenth. 

■n  n  t?  v-6  and  E-8,  some  esti- 
prom  the  results  simphtied  method  cl  corn- 

mate  can  be  made  ol  the  ‘  ^  ^  1„  section  IIA.  In  each 

putlng  pressure  rise  in  a  1  ng  r  ^  ^ 


-  -i  Piire 


X.  Pounds  per  square  inch,  absolute 
t  Pounds  per  square  inch,  gauge. 
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tn  ^  nrcSBurc  equal  to  outside 

ovo  time  and  zero  overpressure  »  iill  pressure  and 

rise  fi'om  zeio  time  overestimates  the  m 

"ossur.  at  a  tl»a  equal  to  V/^  to  «lyo  to  ults 

Lnuorostlmates  the  ""  ^“uon  ant.  ..caauro»c,nt  in  the 

in  good  oeroe»<ent  with  °  '  g  „ithougli  the  prosonco  ol  ” 

cube  model  reported  in  1  igu  pressure  record  m  the  larq 

dilations  (caused  by  It  should  be  noted  also  tha 

„odol  mauea  cleat  blalt  that  is  occupied  by  the  opening 

(rpening  Iraction)  is  over  9  9=-“*  “  ^/g„alght  lino  estimate 

‘rc^nt  in  tho  small  lour-ln  in  yigutc  E-1  ior 

-^\":::ersr;rrs\l“rdrat:Tnroponi„g  traction  is  al„ost  e  percent. 

”  St  that  the  simplified  method  of  esti- 

Tbo  data  Iro»  tho  models  adequate  when 

mating  pressure  rise  in  a  room  sot  tor  ginst  ™vo  As 

outside  pressure  decay  Is  slow  important  diffraction  P^asc  i 

dll  be  noted  later,  wi  structure,  more  sophistical 

the  blast  wave  interaction  with  tn 

may  be  justified.  ; 

2.  outflow.  All  the  “'““'^“^/p'rsturf afoa"'-  “ 
room  pressure  eventually  field  and  if  the  fall  of 

all  openings  are  into  the  same  outsi  P  is  not  likely  to  he 

oitside  pressure  Is  dteahy  and  sl»w  this  ^ 

Of  practical  importance.  outside  pressure  histories,  a 

openings  that  are  affected  by  on  one  or  more  walls  of  the  ro  . 

Significant  outward  O-^dSsure  may  dove  p  npa„ 

TO  calculate  this  pressure,  the  rate 
Ine  must  be  found. 

.  .  the  exterior  atmosphere  and  the  aperture 
TO  compute  inflow  we  treated  ^  ^^ioh  quasi¬ 
leading  into  the  room  as  increase  within  the  room  was  corn- 

steady  flow  was  maintained.  o  „  of  mass  and  energy  from  the  larg 

pr/ed  as  resulting  from  the  treated  in  the  same  way 

outside  reservoir  is  opposite.  For  outflow,  Eqs.  (  )  an 

orcept  that  the  ““the  duet  now  originates  in  the  room  and 

(2)  Still  apply,  c  t  be  replaced  by 

the  adiabatic  expansion  law,  Eq.  (3).  m 

(16) 

^2  "  ^3  P'  J 

^  wnrdrt  to  form  an  equation  identical 
^2)  and  (16)  may  be  combined 
Equations  (1)  •  f  ^ 

to  Eq.  (5)  except  that  now 


f 


B  = 


2Y 

V  +  1 


p' 
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But,  as  before, 


and 


A 


Y  -  1 

Y  +  1 


In  this  form,  Equation  (5)  has  two  solutions  in  the  range  1  >  y  >  0 
whenever  y  == 


1 


or, 


P  P  ' 

1  3 

Pl"  ^  P3 


Y 


Those  two  solutions  merge  into  a  single  solution  at  y  - 


(17) 


if 


(0.9094) 


Furthermore,  if  y  =  1.4,  Equation  (5)  has  two  solutions  but  only  one  in 
the  range  of  0  ^  y  <  1  whenever 


P'  P  P 

-4^^X0.9094)  ^ 


(18) 


Whenever 


P  P 

An  <  (0.9094)  Aw 
D.’  Po  ■ 


(19) 


no  solution  to  Equation  (5)  exists. 


Since  the  quantity  "Y  is  constant  along  an  isentrope  and  increases 
with  entropy,  when  Inequality  (17)  is  true  during  outflow,  specific 
entropy  outside  is  greater  than  or  equal  to  specific  entropy  within  the 

room. 


The  specific  entropy  (i.e.,  entropy  per  unit  mass)  in  the  room  at 
any  time  during  filling  can  be  formally  calculated  as  follows. 

First  we  note  that  temperature  of  air  in  the  room  rises  during  the 
whole  inflow  period,  as  can  be  seen'by  treating  Am  as  a  true  differential 
and  writing  Equation  (12)  as. 
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•  dm 
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anvi  dividing*  numerator  and  denominator  by  V,  and  noting  that  density 
m.  <  d 

,  we  find 

p  Y  T  -  T 

In  other  words  the  process  of  room  filling  during  the  time  period  t  to 

t  +  dt  results  in  a  temperature  increase  dT  that  bears  the  relation 

3 


(dT  )  V  T 

3  isentropic  ^  3 


Y  T  -  T 

^  1  3 

Y  -  T 


to  the  temperature  increase  in  an  isentropic  process  providing  the  same 


density  increase.  When  , 


^^^3) isentropic 


Now  the  entropy  change  dS  in  a  process  at  constant  density  is 


where  c^  is  specific  heat  at  constant  volume.  If  we  imagine  a  return  to 
the  isentrope  (after  an  increment  of  filling)  by  a  process  at  constant 
density  or  volume,  the  specific  entropy  change  within  the  room  resulting 
from  the  filling  increment  becomes 


=  — 
3  T 


dT  -  (dT  ) 

3  3  isentropic 


dS  =  — 
3 


Y  T  - 

V  ^  12 


Y  T  -  T 
^  1  3 


which  on  substitution  of  the  quantity  dT^  fi^om  Eq.  (20)  bocomos 


^v  ^ 

ds  =  — —  (T  -  T,  )dm 
3  T  m  1  3' 

3  3 


Therefore,  during  inflow  (i.e.,  when  dm  >  0)  specific  entropy  within  the 
room  will  increase  until  T3  =  or  whichever  occurs  first. 
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py.  ncnce,  initially  boloro  inllow  begins 


so  created,  Inemalitv  flg^  ‘’°°”'  ‘“'"“f'*’  '"ti-opy  be 

the  outflo;  e^ioL  »i!e  Te  Pt 

cific  entropy  within  the  roL"  °  Tb®  tlsb  in  spe 

py  wixnin  the  room  as  a  result  of  inflow  is: 


p  =P 

/3  1  T  -  T 

■—-ot 

.0  Vs  •“ 


wi^ktew  tharr  °rlsla’  tl«s.  Froi.  Eq.  (20) 

»ay  “Z  thus  ^  ^ 

Tl(t),  tJ)  anr»  rtn  te  ""r  0P«t»l"ss 

Will  be  rwLd  When'  *“'’1”  “e  room 


’'3  “’’i 


If  we  assume  that  at  this  tine  inflow  is  still  under  way,  then 


P,  >  P 
1  3 


which  implies, when  y  >  1,  that 


1-i  1-i 

P  ^  >  P  Y 

1  3 

Y“1  ^-1 


Now  from  the  equality  of  the  temperature 
law  implies 


s  at  this  time,  the  perfect 


!i  !3 

“^3 


Multiplying  both  sides  by  the  last  inequality 


above,  we  find 
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>  Y  p  Y 
1  <-1 


Y  Y 

P  P 

3 

Hence  in  general,  the  possibility  may  not  be  ruled  out  that  Inequality 
(19)  will  be  satisfied.  ^Vhether  the  outflow  equations  as  presented  here 
have  a  solution  will  depend  on  the  nature  of  the  function  P^(t).  If  out¬ 
flow  is  to  be  calculated,  then  to  provide  an  initial  disparity  between 
specific  entropy  outside  and  inside,  the  inflow  period  must  be  treate 
using  the  initial  density  behind  the  shock  front  and  computed  from  the 

Hugoniot  relation: 

(Y+1)P^^+  (Y-l)P^ 

^lo"^o  (y-1)Pj^q+  (Y+DPq 

subsequent  outside  air  densities  may  be  calculated  from  the  adiabatic 
law:  •  ^ 


in  Eqs.  (23)  and  (24),  and  P^^  refer  to  the  air  density  and  ab¬ 

solute  pressure  immediately  behind  the  shock  front.  When  peak  pressure 
at  an  opening  is  reached  by  reflection  of  an  incident  shock  wave  from  a 
wall  Eq.  (24)  is  not  correct  regardless  of  whether  P^^  is  taken  as  peak 
incident  absolute  pressure  or  peak  reflected  absolute  pressure,  but  the 
error  in  using  Eq.  (24)  is  small.  We  will  arbitrarily  consider  and 

P,  as  representing  conditions  behind  the  free-field  shock  front.  Use 
of°Eqs.  (23)  and  (24)  may  make  it  possible  to  satisfy  the  reverse  of 

Inequality  (19). 

Should  such  still  not  be  possible,  other  calculational  methods^must 
be  used  for  the  outflow  phase,  such  as  that  proposed  by  Melichar  »  or 
that  reported  by  IIT  Research  Institute. These  methods  equate  duct 
pressure  Po  with  outside  pressure  during  outflow;  the  IIT  investigators 
then  fit  observed  outflow  pressure  data  by  choosing  values  for  the  is- 
charge  coefficient  and  for  the  ratio  of  inside  to  outside  pressure  at 

the  time  of  flow  reversal. 
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Hand  Calculation.  In  constructing  from  the  foregoing 
equations  a  calculational  scheme  for  estimating  the  parameters  of  flow 
into  and  out  of  a  single  room  with  a  single  opening  wo  start  with  a  scries 
of  values  of  P^,  one  for  each  time  step.  These  may  bo  obtained  by  linear 
interpolation  from  a  given  table  of  outside  pressure  as  a  function  of  time 
and  each  value  should  pertain  to  the  center  of  the  time  interval.  The 
size  of  the  time  interval,  At,  itself  is  arbitrary,  but  it  should  be  no 
greater  than  the  quantity  T;  presumably  up  to  a  limit,  greater  accuracy 
results  from  smaller  values  of  At.  The  size  of  At  may  bo  changed  during 
the  calculation  when  the  rate  of  change  of  flow  parameters  changes.  Wo 
also  need  values  for  ambient  pressui'e  and  density  p 

Ttt'o  methods  of  calculation  are  shown  below.  The  first  is  that  used 
to  produce  Curves  F  in  Figures  E-4,  E-5,  E-6,  and  E-8,  namely,  that  based 
on  Eqs.  (1),  (2)  and  (3)  for  inflow  or  Eqs.  (1),  (2)  and  (24)  for  outflow, 
and  in  the  outline  below  it  is  called  Method  F.  This  method  has  the  ad¬ 
vantage  of  great  simplicity  and  of  not  requiring  knowledge  of  empirical 
constants;  however,  as  will  be  explained  later,  values  of  wind  speed  and 
dynamic  pressure  computed  by  it  are  subject  to  doubt  in  some  cases  and 
for  that  reason  a  method  given  by  IIT  Research  Institute  is  included  also. 
The  latter  method  is  responsible  for  Curves  D  in  the  figures  and,  there¬ 
fore,  m  the  outline  below  it  is  called  Method  D.  As  noted  earlier 
Method  D  for  the  unchoked  flow  case  is  numerically  equivalent  to  the 
calculation  used  by  Melichar.^^i 

Average  pressure  inside  a  room  and  dynamic  pressure  in  the  single 
opening  to  an  outside  reservoir  whose  pressure  variation  in  time  is  known 
may  be  calculated  as  functions  of  time  by  the  sequential  application  of 
the  steps  stated  below.  Each  cycle  through  a  series  of  steps  completes 
the  calculation  for  one  time  interval.  The  first  three  steps  are  executed 
only  during  the  first  cycle;  subsequent  passes  begin  with  step  (4),  as  in¬ 
dicated  in  the  outline.  Step  (5)  is  a  branch  point  to  separate  sequences 
for  inflow  and  outflow,  chosen  according  to  a  criterion  given  in  step  (5). 
There  are  further  branches  (a)  to  Method  D  or  F,  chosen  at  the  discretion 
of  the  user  at  each  time,  and  (b)  under  Method  D  to  choked  or  unchoked 
flow,  determined  by  stated  criteria.  Throughout  the  outline,  the  quantity 
Y  has  been  set  equal  to  1.4.  ■ 


(1) 

Set  P'  =  P  and 

3  o 

p;  = 

O  o 

and  t  =  0. 

(2) 

Compute 

n  —  n 

6P^  +  P 

lo  o 

lo  .  o 

P,  +  6P 
lo  ■  o 

where  P  is  the 

shock  front* 

absolute 

pressure  immediately  behind  the 
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(3)  Choose  value  of  At  (see  opening  paragraph  of  this  Section  IIB3). 

(4)  Determine  outside  pressure  at  the  current  time;  i.e.,  deter¬ 
mine  from  the  known  reservoir  pressure  history  (pressure 

variation  with  time).  During  first  time  interval 

(5)  Determine  direction  of  flow;  i.e.,  if  >  P3,  flow  is  inward; 
go  to  step  (6).  Otherwise  flow  is  outward;  go  to  step  (28). 


Inflow 

(6)  Compute 

Branch  to  selected  Method  below  for  step  number  (7D)  (Method  D) 
or  step  (7F)  (Method  F).. 

Method  D  (Inflow) 

(7D)  If  ^  0.5283  inflow  is  choked;  go  to  step  (8D).  Otherwise 

inflow  is  unchoked;  go  to  step  (18D). 


Choked  Inflow 


2.4 


(8D)  Am 


choked 


1/2 


/  2  \0.4  ^ 


r  -lV2 

=  0.6847  K  [p  P  f]  a  At 
X  1  ^ 


At 

2 


Using  the  recommended  value^  of  K  =  0,70  this  becomes 

1/2 

^  1  ^  =  0.4793  [  f]  At 

choked  11  2 


(9D)  P2=P; 


(lOD) 


,0.7143 


IJ 


*  The  factor  F  will  often  be  necessary  for  consistency  of  units.  For  ^ 
example,  if  U2  is  to  be  in  ft/sec  and  P^^  is  in  Ib/in^  and  pj^  in  Ib/ft  , 
then  F  =  32. X  144.  =  4,608. 
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(25D)  P'  =  P3 
(26D)  P'  =  P3 

(27D)  Advance  time  by  amount  At  and  return  to  step  (4) . 
Method  F  (Inflow) 


(7F) 

=  0.1912  P, 

2  1 

(8F) 

^2  =  Pi 

—  =  0.3067  p. 

(9F) 

2  1 

^  =  ^  1 

r!i.!2l.2.637r^lF 

Lpi  P2J  ["iJ 

(lOF) 

1/2 

p^A^  At  =  0.498  (P^P^F)  A^At 

(IIF) 

^3=^3 

P  Am 

1  3 

+  1.4  — 

P  V 

^1  3 

Am 

0 

(12F) 

P3  =  P3 

0 

(13F) 

If  desired,  the  dynamic  pressure,  q^,  in 

the  doorway  may 

be  calculated; 

(14F) 

1 

^2*2* 

2 

^2 

3  -  =  0.4044  P, 

2  F  1 

(15F) 

P3=  P3 

(16F) 

^3 

(17F) 

Advance 

time  by  amount  At  and  return  to 

step  (4) . 

Outflow 

(28) 

Branch 

to  selected  Method  below  for  step 

(29D)  (Method  D) 

or  step 

(29F)  (Method  F). 
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Me thod  D  (Outflow) 


^29D')  If  Pi/P2  ^  0.5283,  outflow  is  choked;  go  to  step  (SOD), 
Otherwise  outflow  is  unchokod;  go  to  step  (39D) . 


Choked  Outflow 


1/2 

'^Cl.o.od  = 


For  K  ==  1.0  this  becomes 


17 


1/2 

=  -0.6847  [p'  P'  f3  A^At 
chokod  ^  ^ 


(31D)  Pg  = 


(32D)  Po  =  P3 


0.7143 


^choked 

"2-itA2  PjK 


(34D)  Pr,  = 


^Choked  ^3 


Am 


(35D)  P3  =  P3  ■*■  7 


•  (36D)  P3  =  P3 

(37D)  P'  =  P3 

(38D)  Advance  time  by ' amount  At  and  return  to  step  (4) 


Tlnchoked  Outflow 


(39D)  P„  =  P. 


(40D)  P2  =  P' 


(41D)  u  =  7 


0.7143 


P' 

L  3  J 


!i.!2 

"1  1=2 
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'^"unohoted  “  -‘’2  “2  \ 


(43D)  ^3=^3+ 


Am  u  1  ^ 

unchoked  d 

V  p' 

3  3 


Am 


(44D)  P3  =  P3  + 


unchoked 


(45D)  P'  =  Pg 
(46D)  p'  =  pg 


(47D)  Advance  time  by  amount  At  and  return  to  step  (4). 


Method  F  (Outflow) 


7  ^3 

(29F)  B  =  -  — 

^  /  6  p 


/,  0.7143 


L^J 


(30F)  Solve 


0.7143  ,  . 

By  =  y  +  A  for  y. 


(31F)  P^  =  yPj^ 


(32F)  P^  =  P' 


.  P' 

L  3 


0.7143 


(33F)  ^  =  -P2  “2  \ 


1  Am 

(34F)  P3  =  P'  +  1.4  -  - 

1  3 


Am 


(35F  P3  =  P3  V 


(36F)  P'  =  P3 


(37r)  p;  =  P3 


(38F)  Advance  time  by  amount  At  and  return  to  step  (4) . 
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.  1  diflorencos  WWoen  the  ““ 

thooretxca^^  *'“'%ror Juc.,^.oht  between  the.. 

„e  stribing.  permit  “  •  ^j,„e  is  lone  comperec 

dence  from  slowly,  then  the  analogy  under 

11  the  room  soonO  acroos  the  X  about  the  cor- 

,1th  "“"“should  be  good;  but  even  so  th^  "‘  r. 

lyibS  «■=""“”*  discharge  coeflicient,  K  the  diller 

rect  value  o  pressure  is  large,  filling  is  rapid,  P 

on  the  oalculat  ^  Figure  E-4.  seen  in  the  lol- 

enca  beween  0“"®“  preterred.  However,  as  w  ^^p<,  speed 

haps  Method  r  should  ^  pr  ^p  s„pmalous 

lowing  paragraphs,  this 
through  the  opening. 

n  Wind  Speed.-:i.SHS^ -  ,  p  ppos- 

C.  —  A  dvnamic  prcb 

eters  of  interest  are  wind  parameters  only 

rpn^b"-  -  .  p  ^  ppp  son  dFnamio  pressures  in 

„„f ortunately , 

shoch-lillidB  "  rhe  predictions  ol  ““  ®p  pressure  are  otten 

■‘^‘r::rdsTrriarge:  but  the  vajues  PP^PP,S  IP  the 

rnTaragreement.  -  -rirthfoSy  pertinent  parameter. 

Stream,  the  dynamic  p  calculated  from 

outflow  wind  speed  V  Dynamic 

During  Jentropic  assumption  is  made.  q. 

Eq.  (9)  "'^\^3tomarily  defined  as 

pressure,  Q2’ 


u. 
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Using  Eq.  (9)  this  can  be  written: 

^2  ^  I  \  ■  (25) 

Because  of  the  neglect  of  transient  effects,  Eqs .  (9)  and  (25)  do  not  re¬ 
flect  the  influence  of  changes  in  the  room  pressure,  and  values  calculated 
from  them  may  be  unrealistic.  Peak  speeds  found  from  Eq.  (9)  for  the  ex¬ 
periments  reported  in  Figures  E-4,  E-5,  E-6,  and  E-8  are  supersonic,  i,e., 
above  local  sound  speed  in  the  doorway. 

Because  of  its  use  of  the  choking  analogy  Method  D  never  calculates 
supersonic  speeds;  but  it  does  neglect  transient  shock  fronts  which  could 
alter  the  flow  substantially  in  fast  filling  rooms. 

Table  E-1  contains  the  several  calculated  values  of  peak  speed  and 
dynamic  pressure  corresponding  to  each  of  the  experiments  reported  in 
Figures  E-4 ,  E-5,  E-6,  and  E-8.  The  table  suggests  that  use  of  Method  F 
may  result  in  overestimate  of  that  important  criterion  for  damage  potential, 
dynamic  pressure.  As  noted  above,  because  signals  from  the  filling  room 
into  the  outside  reservoir  are  neglected,  Method  F  in  addition  to  over¬ 
estimating  peak  dynamic  pressure  may  overestimate  the  time-average  dynamic 
pressure  in  the  doorway. 

There  are  no  direct  observations  of  wind  speed  in  the  opening  to  com¬ 
pare  with  calculations.  Coulter,®  however,  had  reported  (1)  measurements 
of  the  acceleration  of  an  g-inch  diameter  nylon  ball  placed  in  the  doorway 
of  a  small  model  struck  head-on  by  a  weak  shock  wave,  and  (2)  measurements 
of  the  motion  of  smoke  columns  inside  these  models  during  shock  filling. 

With  the  help  of  acceleration  coefficients  derived  (theoretically  and  ex¬ 
perimentally)  by  Bowen  and  others^^  we  can  compute  motion  of  the  nylon 
sphere  from  Eq .  (25).  Coulter  allowed  a  shock  front  of  overpressure  equal 
to  4.89  psig  to  strike  a  reflecting  plate  in  which  an  entrance  1x4  inches 
was  cut.  Since  the  chamber  behind  the  opening  had  a  volume  to  opening  area 
ratio  (V/A)  equal  to  1.33,  filling  was  essentially  complete  in  1.33/2  = 
0.667  ms,  during  which  time  there  was  little  or  no  decay  of  the  incident 
wave.  If  we  assume  ambient  pressure  equal  to  14.7  psi,  reflection  at 
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Table  E-1 

PEAK  VALUES  OF  OTHER  FLOW  PARAMETERS 


C'X- 

F 

G 

Experiment 

S 

^^2 

Fig.  E-4 

7.98 

1045. 

9.93 

1650. 

8.02 

985. 

Fig.  E-5 

10.2 

1180. , 

11.2 

1680. 

10.2 

1110. 

Fig.  E-6 

7.79 

972. 

9.74 

1640. 

7.80 

972. 

Fig.  E-8 

4.08 

747. 

7.37 

1634. 

*  Correct  adiabatic  law  has  been  substituted  for  the  inverted  form 
appearing  in  Melichar's  FORAST  code.^^ 


normal  incidence  of  a  4.89  psi  shock  produces  a  reflected  overpressure 
equal  to  11.1  psi.^  Hence  the  value  of  P^  in  Eq .  (25)  becomes  11.1  +  14.7 
25.8  psi,  and 

q  =  -  X  25.8  (1  -  0.1912) 

2  2 

=  10.5  psi 

2 

or  q  =  48100  Ib/sec  ft 

2 

In  this  expression  "lb"  is  a  unit  of  mass. 

A  still  object  in  a  stream  of  moving  air  is  accelerated  according  to 
the  formula 


where  a  is  an  acceleration  coefficient  characteristic  of  the  object  and 
q  is  the  dynamic  pressure  as  defined  here.  Eq.  (26)  continues  to  hold 
for  the  moving  object  as  long  as  its  speed  is  small  compared  to  air  speed. 
Bowen  and  other^^  suggest  a  value  of  o'  =  0.138  ft^/lb  for  a  1/8  inch 
diameter  steel  sphere.  Since  a  is  inversely  proportional  to  mass,  the 
corresponding  value  for  a  nylon  sphere  of  the  same  size  is  approximately 
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2 

7.87  X  0,138  =  1.09  ft  /lb  (specific  gravity  of  iron  equals  7,87).  Hence 
nv'v'el<'rntlon  of  the  nylon  sphere  should  be,  according  to  Eq .  (26), 

—  =  1.09  X  48100  =  52200  ft/soc^ 

CJ  u 

and  the  final  speed  of  the  sphere  when  filling  is  complete  667  p,sec  after 
impact  becomes : 


V  =  52200  X  667  x  lo"®  =  34.8  ft/sec 

1  2  —12 

and  the  displacement  is  s  =  X  52200  X  (667)  X  10  =  0.0116  ft  = 

0.14  inch.  Coulter  measured  terminal  speed  after  700  y,sec  as  29,8  ft/sec. 
Observed  displacement  in  this  same  time  was  approximately  0.10  inch,  which 
means  the  ball  did  not  move  far  from  the  center  of  the  doorway  during  the 
period  of  observation  and  was  subjected  to  only  slightly  attenuated  peak 
dynamic  pressure . 

Corresponding  to  Eq.  (25),  Method  D  produces  the  following  equation 
for  the  calculation  of  dynamic  pressure  in  the  opening: 

1 


The  value  of  q^  in  Eq,  (27)  is  not  constant  in  time  even  when  does  not 
change  appreciably.  For  simplicity,  we  may  assume  P3'  varies  linearly  from 
14.7  psia  to  P^  over  the  filling  time  0.667  ms.  Again,  we  take  P^^  =  25.8 
psia  and  find  the  following  variation  of  dynamic  pressure  with  time; 


time  (ms)  0 

0.2 

0.4 

0.6 

0.667 

Pg'  (psia)  14.7 

18.0 

21.35 

24.7 

25.8 

q^  (psi)  8.97 

6.80 

4.15 

1.08 

0.00 

The  tabulation  shows  a  nearly  linear  variation  of  q^.  Integrating  Eq.  (26) 
approximately,  we  find  the  final  value  of  speed  to  be: 
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-  X  8.97  —2  ^ 

lb  2 


.667  X  10"^  soc  144  32 

i 


=  15.0 


n  Tv,«  rplativelv  better  agreement 

:«reK;eri:rnt''fuLrsLrr;;:t;c:  p  ;han  ^  -y  r.nect  t.o 

fact  that  the  shpok  lilime  In  this  case  is  last. 

However,  it  would  seem  that  any  wariortultous,  since 

oalculated  and  “ch  ^ 

Coulter’s  photographs  of  events  calculated  filling 

dilfraotlon  episode  very  near  y  co  observation  was  more  a  shock 

rime,  that  is,  air  flow  durin.  Some  support 

related  process  than  a  case  q  observations  reported  by 

for  this  View  comes  fro.  the  -«nd  „,od  for  measure- 

Coulter:  motion  of  smoke  in  Coulter’s  observed  values 

ments  of  the  motion  of  the  Lion  are  one  or  two 

of  wind  speed  and  dynamic  pressu  values  calculated  by  the  methods 

orders  of  magnitude  lower  than  compara  nethod  measure 

presented  here.  However  “  on  his  observations 

rrfLstrSLrLtL  Lre:-a::  Lther  from  the  core  of  the  flow. 
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To  provide  a  useful  appraisal  of  the  effects  of  inflowing  air  within 
an  open  shelter,  the  attenuatxon  of  wind  speed  and  the  snreadin.^^  of  the 
•  ■ciS  W..G’  » r. T'w — ^ tr,c-  ciGc-yc-r*  ■t’r.C'  Gr-clcsurc  '-"■'♦c** 

be  accounted  for.  As  in  our  description  of  the  averag:e  pressure  rise 
within  the  room,  reliance  will  again  bo  placed  on  a  Quasi  “-static  approach 
with  account  taken  of  the  transient  episode  before  the  establishment  of 
adiabatic  flow  by  means  of  a  semiempirical  correction. 

The  steady  jet  described  by  Abramovich^^  is  illustrated  in  Figure 
where  air  is  considered  to  be  flowing  from  a  reservoir  at  high  pressure  on 
the  left  through  the  aperture  into  the  region  of  low  pressure  (i.e.,  the 
interior  of  the  shelter)  on  the  right.  In  the  theory,  the  entry  aperture 
may  be  a  slit  of  essentially  infinite  length  and  width  or  it  may  be  a 
circular  opening  of  diameter  b^.  Although  the  description  of  the  jet  is 
different  for  the  slit  than  for  the  circle,  the  difference  is  not  impor¬ 
tant  here.  Since  most  actual  openings  will  be  somewhere  between  the  two, 
it  will  be  assumed  that  any  actual  opening  has  been  approximated  by  a 
slit  or  by  a  circle  and  only  one  set  of  formulas  will  be  given.  As  shown 
in  Figure  E-9A,  inflowing  air  generally  fans  out  into  the  receiving  res¬ 
ervoir,  slowing  down  after  passing  through  a  conical  or  wedge-shaped 
”core,*^  Air  speed  is  constant  everywhere  within  the  core. 


Abramovich  characterizes  the  jet  by  three  parameters:  dimension 
of  the  opening  bo  ;  air  speed  in  the  axial  direction  at  the  initial  cross 
section  Uq  ;  and  ratio  9  between  temperature  in  the  core  and  that  in  the 
receiving  reservoir  (i.e.,  the  shelter  space).  The  flow  is  driven  by  a 
pressure  difference  between  the  two  reservoirs,  but  it  appears  to  be  a 
fact  that  pressure  within  the  usual  subsonic  jet  itself  is  uniform  and 
equal  to  that  in  the  receiving  reservoir.  In  other  words,  pressure 
equilibrium  is  ,quickly  established  between  the  jet  and  the  air  already 
in  the  reservoir.  In  this  case,  where  the  air  in  the  high  pressure  res¬ 
ervoir  has  been  shock  compressed  from  a  state  identical  to  that  origi¬ 
nally  in  the  low  pressure  reservoir,  the  existence  of  pressure  equilib¬ 
rium  between  the  inflowing  air  and  interior  air  requires  that  density 
and  temperature  in  the  jet  nearly  equal  density  and  temperature  in  the 
room.  The  slight  difference  is  attributable  to  the  fact  that  outside 
air  has  passed  through  the  shock  front  and  acquired  entropy  compared 
with  the  unshocked  or  only  weakly  shocked  air  within  the  room.  Such 
difference  is  entirely  negligible  for  present  purposes.  Thus,  in  this 
discussion,  the  parameter  0  equals  1. 

Furthermore,  if  pressure  is  uniform  within  the  jet,  the  flux  of 
axial  momentum  must  be  the  same  at  each  cross  section  of  the  jet,  that 
is : 
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FIGURE  E-9A  SCHEMATIC  OF  JET  FLOW 
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By  making  use  of  a  control  surface  passing  through  the  aperture  (shown 
by  the  dashed  line  in  Figure  E-9B) ,  the  value  of  Uq  in  Equation  28  can  be 
calculated,  since  by  conservation  of  momentum 


flux  =  (P  -  P  ) A 
1  3  2 

Thus, 


(29) 


Alternatively,  by  Equation  2 


P  -  P 
1  2 


or 


P.  = 


'‘2''2  =  ^“2 


"“2  ^  ‘’2 


which  can  be  substituted  for  P-j_  in  Equation  29  to  yield: 
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This  may  be  a  useful  lower  bound  since  P^  is  generally  small.  Both 
Equations  29  and  30  show  that  u^  falls  as  the  room  fills. 
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Expcriiricntal  observations'^^  in  model  rooms  and  detailed  nurnerJcal 
simulations'^^' of  shock  room  filling  both  dornonstrato  that  in  room 
filling  streams  dynamic  pressure  roaches  values  about  tv/ice  those  given 
by  Equation  31,  Thus^  instead  of  Kquations  nitd  'At  t}//v  bs 

recommended  on  empirical  grounds : 


u 
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(32) 
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cox*e 


(33) 


In  most  cases  the  relation 


q 


core 


(34) 


has  been  found  empirically  to  be  approximately  true  at  the  beginning  of 
the  jet  period.  The  use  of  Equation  34  throughout  the  filling  period 
appears  from  experimental  evidence  to  be  highly  conservative,  since  in 
addition  to  the  spatial  attenuation,  observations  reveal  a  fall  in  core 
dynamic  pressure  with  time  during  the  equalization  of  inside  and  outside 
pressures . 

The  reason  for  the  failure  of  Equations  29  and  31  is  not  understood; 
however,  with  further  empirical  corrections,  steady-state  jet  theory 
may  be  applied  with  fair  success  to  describe  the  room  filling  airstream 
quantitatively. 

The  first  departure  from  straightforward  steady  jet  behavior  is 
observed  during  the  period  immediately  after  shock  arrival  at  the  open¬ 
ing,  when  air  motion  in  the  neighborhood  of  the  opening  is  governed  by 
the  laws  of  shock  waves.  At  that  time,  particle  speed  is  relatively  slow 
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but  as  the  diffracted  waves  reverberate  in  the  opening,  conditions  for 
jet-like  flow  are  established.  This  build-up  of  flow  through  the  open¬ 
ing  lasts  a  length  of  time 


(35) 


where  c  is  the  speed  of  sound.  If  this  time  interval  is  long  enough  to 
be  important,  it  is  recommended  that  dynamic  pressure  in  the  opening  be 
treated  as  varying  linearly  in  time  during  to  from  cjq  =  0  at  shock  arrival 
(t  =  O)  to  qo  =  Pi  -  P3  at  the  end  of  the  diffraction  phase  (t  =  to). 
During  the  time  t  <  to  ,  the  variation  in  P3  with  time  cannot  be  calcu¬ 
lated  since  the  filling  law  during  diffraction  at  the  opening  is  extremely 
complicated.  If  substantial  filling  is  thought  to  take  place  during  the 
period  t  <  t^  ,  the  inflow  rate  at  any  time  step  may  be  found  from  the 
formula: 


Am„ 
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P  u  A„ 
^3  o  2 


=  A, 


o  3 


1/2 


(36) 


The  jet  is  not  established  simultaneously  throughout  its  final 
length.  An  estimate  of  the  time  between  shock  arrival  and  the  appear 
ance  of  jet-like  flow  at  an  interior  point  i  appears  to  be: 


(37) 


where  is  the  distance  from  the  opening  to  the  point.  Again  a  fair 
appraisax  of  the  time  variation  of  dynamic  pressure  at  point  i  during 
this  preliminary  period  would  seem  to  be  a  straight  line,  beginning  at 
the  moment  of  shock  arrival  at  point  i . 

If  the  point  i  is  near  a  wall  blocking  flow  from  the  opening,  then 
dynamic  pressure  there  may  never  reach  a  value  expected  in  a  fully  devel 
oped  flow  at  that  distance  from  the  opening. 


The  fully  developed  steady  subsonic  ict  ,t 

ins  throe  sets  of  formulas,  each  sot  valid  la  one  ho  t hror*’ 
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.  the  transition  region,  there  is  no  core; 
..  ».uc.  a.  .al.  an.^a  eaaaX  ao 


a  Si  12°30 ' 


or  slope 


-  =  0.22 

X 


where  b  is  now 


the  half  width  of  the  jet. 

is  at  a  maximum  at  the  jet  axis  and 
In  any  cross  section,  air  between  at  a  perpendicular  dis- 


in  any  -  -  In  between 

falls  to  zero  at  the  outer  boundary. 

tance  y  <  b  from  the  axis 


1  5  2 

„  =  u  (1  - 

m 


where  u„  =  speed  along 


^  n  _  I  as  before,  except  that  y  and  b 
the  aad  T|  -  ^ 


have  alightly  dlfierent  Beahinge,  as  shovm  in 
Sketched  in  Figure  E-9A. 
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Any  potential  for  damage 
from  Equation  29.  Changes  in 
in  that  equation  are  not  felt 


ding  in  the  airstream  stems  ultimately 
values  of  the  thermodynamic 
ndiately  and  simultaneously  throug 
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the  strofim.  However,  the  assumption  of  instantaneous  response  will  fall 
to  be  conservative  only  after  the  end  of  inflow,  i.c,,  wlicn  and 

Uo  as  calculated  by  Equation  29  becomes  0.  In  fact,  the  stream  cannot 
.stop  abruptly.  If  the  airstroam  deep  within  a  shelter  remains  dangerous 
as  filling  concludes,  then  a  linear  decay  in  dynamic  pressure  at  the  in¬ 
terior  location  may  be  assumed,  lasting 


or,  more  simply,  inertial  effects  during  decay  may  be  accounted  for  by 
neglecting  inertial  effects  during  build-up  of  the  stream  at  the  interior 
point . 

Obstacles  to  flow,  such  as  corners  or  barriers,  may  be  found  in  a 
shelter  or  may  be  deliberately  designed  into  the  structure.  As  far  as 
the  main  shelter  area  is  concerned,  such  entry  barriers  or  mazes  have 
the  effect  of  increasing  the  length  of  the  diffraction  phase.  Thus  if 
the  driving  pressure  is  falling  rapidly  enough,  the  duration  and/or  in¬ 
tensity  of  the  airstream  flowing  into  the  shelter  may  be  reduced  by  the 
presence  of  a  simple  barrier  at  the  entry,  as  sketched  in  Figure  E-9C. 
Increase  in  the  duration  of  the  diffraction  phase  (i.e.,  the  time  inter¬ 
val  t^)  is  a  possible  explanation  for  the  positive  results  of  Coulter's 
barrier  experiments.^’’  However,  Equations  32  and  33  appear  to  offer  the 
possibility  of  baffles  that  would  with  a  sacrifice  of  shelter  space  greatly 
reduce  the  hazard  connected  with  the  filling  sti’oam.  For  example,  a  hold¬ 
ing  i-oom  or  foyer  is  shown  in  Figure  E-9D  that  could  conceivably  bo  de¬ 
signed  so  that  prossuro  within  it  during  most  of  tho  fill  period  would 
bo  loss  tlian  that  present  outside  tlio  slioltor.  If  sucli  a  design  is  pos¬ 
sible,  Equation  33  tolls  us  that  the  dynamic  prossviro  in  tho  core  of  the 
jot  emerging  into  tho  actual  shelter  space  would  bo  loss  than  it  would 
bo  without  tho  foyer.  The  difference  between  the  holding  rooms  shown  in 
Figures  E-9D  and  E-9E  should  be  noted.  The  driving  pressure  for  the  jet 
flowing  into  the  shelter  space  of  Figure  E-9E  is  essentially  the  outside 
pressure,  not  that  in  the  holding  room.  This  arrangement  benefits'  only 
from  the  increased  delay  tQ  .  The  hazardous  jet  of  Figure  E-9D  on  the 
other  hand  is  driven  by  the  pressure  in  the  holding  room.  If  the  foyer 
volume  V  and  door  areas  and  A2  could  be  designed  to  keep  the  average 
pressure  in  V  substantially  below  ,  then  the  arrangement  would  lead 
to  reduced  values  of  Qq  after  the  delay  t^  has  elapsed. 

Some  of  the  same  considerations  may  be  used  in  understanding  stream 
deflection  within  the  shelter  space,  such  as  the  situation  illustrated 
in  Figure  E-9F,  where  the  airstream  flows  into  a  corner  of  the  shelter 
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FIGURE  E-9D  POSSIBLE  HOLDING  ROOM  TO  REDUCE  WIND  HAZARD 
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FIGURE  E-9E  IMPROPERLY  DESIGNED  HOLDING  ROOM 


CD  ^ 


FIGURE  E-9F  DEFLECTION  OF  AIRSTREAM  IN  CORNER 


and  eventually  flows  out  of  the  corner  in  a  direction  perpendicular  to 
the  direction  of  inflow.  After  inflow  into  the  corner  develops,  there 
is  a  time  delay  during  which  static  pressure  builds  up  tliroughout  an 
increasing  volume  of  stagnant  air  V.  This  increasing  static  pressure 
drives  the  outflow  in  ever-increasing  magnitude.  If  inflow  remains 
steady,  presumably  an  equilibrium  is  achieved  such  that  inflowing  mass 
matches  the  mass  flowing  out. 

If  temperature  equilibrium  is  also  reached,  then  kinetic  energy 
of  inflow  must  equal  the  energy  of  outflow.  Thus,  the  net  result  after 
enough  time  has  passed  is  a  pure  deflection  of  the  incoming  stream.  In 
small  shelters  subjected  to  megaton  nuclear  weapon  effects,  this  is  the 
condition  to  be  expected. 

Behavior  of  an  airstream  in  the  corner  illustrated  in  Figure  E-9F 
differs  from  that  suggested  by  Figure  E-9D  in  that  the  areas  A]_  and  A2 
need  not  be  equal,  and  temperature  equilibrium  between  the  holding  room 
and  shelter  space  need  not  be  reached  while  blast  pressures  against  the 
structure  exist. 


III.  Multiple  Rooms 


Although  the  algebra  becomes  increasingly  cumbersome  as  rooms  are 
added,  the  foregoing  principles  of  calculation  can  be  applied  to  a  series 
of  connected  rooms.  For  example,  two  rooms  are  represented  in  Figure  E-11 
for  which  two  control  surfaces  may  be  employed.  Through  the  smaller  the 
momentum  flux  is  due  to  flow  through  the  first  doorway  and  the  resulting 
momentum  equation  may  be  written  as: 


A  (P  -  P  ) 
2  '  1  2 


2 

p  u  A 
^2  2  2 


(32) 
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Using  the  larger  control  surface  wo  find: 
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u"  A 
4  4 


Mv'thovl  oiuploys  ihu  two  o<]Ualitie.s  below  (Juring  filling: 


(33) 


(34) 
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Equations  governing  flow 
can  be  written  as  follows: 


into  the  first  of  the  two  connected  rooms 


P  P 

J-  -  JiL.  i  2 

Y-l  “  Y-1  Pg  2  ^2 


(36) 


(37) 


To  these  two  equations  one  of  the  two  possible  pressui’e  relations 
represented  by  Eqs.  (32),  or  (34)  must  be  added  to  form  a  complete  sys¬ 
tem  from  which  ,  u^ ,  and  for  the  time  interval  under  consideration 
may  be  calculated. 


Under  Methods  F  and  G,  the  equations  governing  flow  into  the  first 
of  two  rooms  are  identical  to  those  from  which  the  filling  of  a  single 
room  is  computed.  Calculation  of  flow  into  the  second  or  inner  room  is 
more  complicated  however. 


Because  observed  pressure  buildup  in  two  connected  rooms  has  not 
generally  followed  simple  calculations  based  on  an  assumption  of  uniform 
pressures  over  the  control  surfaces  of  Figure  E-10,  provision  will  be 
made  in  the  analysis  at  the  outset  for  a  correction  term  A.  The  quan¬ 
tity  A  will  be  chosen  to  reconcile  momentum  balances  based  on  inner  and 
outer  control  surfaces  such  as  those  shown  in  Figures  E-3  and  E-7,  that 
is, 
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FIGURE  E-10  CONTROL  SURFACE  USED  IN  CALCULATING  FLOW  INTO 
SECOND  ROOM 
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•^3  ' 

A  =  (P^  -  P3)  A2 


.ill  take  the  correction  terms,  A,  equal  to  the  pressure 
ri^re^inruracros:  the  duct  times  the  area  ol  the  duct.  Thus,  the  cor¬ 
rected  Eq.  (33)  becomes; 


("3  -  V  "4  “4  ^ 


During  filling  of  the  second  room,  there  the  speed 

moving  the  ^“"°^^^i„^^rarea  of  the  first  room,  then, 
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conservation  of  mass  implies 


P3  ^3  ^3  "  P4  %  % 


and  the  assumption  of  isentropy  leads  to 
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Eqs.  (39),  (40),  (41),  and  (42)  make  it  poss 
the  four  unknowns; 


ible  to  compute  by  Method  F 
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The  resulting  solution  for  P  o.n  k 
in  computing  p  4  pnt  in 


3  Torm  similar  to  that  uj 
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and  B  = 
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then  Eq.  (43)  can  be  wril 


“  Y  =  y  +  A 


where  in  general  A  -  A^'t,^  . 

the  dependence  of  .  on  ;  in  “"ri .0““'  “a  \  «  A^, 


of  Eq.  (45)  13  equivalent  to  nevlectinv  tu 
*.rs^  room  in  Eq.  (45)  _  ^  ‘he  averaee  wind  speed  in  the 

After  the  duct  parameters  i  e 

U3,  have  been  found  from  the  f oregoin^^lo.,!? •  ‘'2’  ^4'  ""4'  IT  desired, 

-tee  P^  and  P^  and  densities 
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where 
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{  J1 - 1  +  i  =  —  (P  -  P') 
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Amc  =  u  p  A  .  At  mass  increment  in  second  room 
4  4  4 
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where 
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Outflow  from  both  the  first  and  second  rooms  is  treated  with  similar 
basic  equations  except  Eq,  (37)  is  replaced  by 


1/Y 


p  =s  p  ' 

^2  ^3 


(50) 


and  Eq.  (42)  by 


P .  \  l/Y 


(51) 


Again,  in  the  computation  of  outflow  parameters  in  the  first  doorway 
(i.e.,  p^,  and  u^)  Eqs.  (32),  (36),  and  (37)  are  used  as  set  forth 
above  when  only  a  single  room  was  considered*  For  outflow  through  the 
second  doorway,  two  cases  should  be  distinguished:  (1)  there  is  inflow 
through  the  first  doorway  and  (2)  flow  through  the  first  doorway  is  out¬ 
ward.  These  conditions  affect  the  form  of  Eqs.  (33)  and  (40).  In  case 
(1),  u^  =  0,  and  Eq.  (40)  simplifies  to: 
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but  Eq.  (33)  must  contain  two  correction  terms,  A^  and 
excess  force  against  the  righthand  and  lefthand  walls, 
the  first  room; 


A  ,  to  account  for 
5 

respectively,  of 
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where  x  and  B  are  given  by  Eqs.  (55)  and  (56),  but  now 
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Substituting  Eq.  (59)  into  Eq.  (61), 


^  Itl  I  (1  _  a^)  +  c?  y| 


A  a  -ind  V  <  1  «  1.  from  which  we  conclude; 

Since,  ordinarily,  »  A^,  and  y  ,  . 


A  -in  raqe  (2)  can  be  found  by  writing 
An  alternative  expression  for  A  m  case 

An  alternati  p  aeainst  which  the  outflow  from  the  second 
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for  which  small  Of  becomes  r* 
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One  final  case  remains  to  be  considered,  namely,  outflow  through  the 
first  doorway  combined  with  inflow  through  the  second.  Conditions  in  the 
second  doorwiy  are  computed  by  solving  Eq..(40)  with  U3  =  0,  -Eq.  (42)  and 

(P^-P')A2+(P3-P4)A4=  ^ 
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These  equations  lead  to  Eq.  (44)  in  which 


P 

4 


Y+1 


1  In  summary  then  we  see  that  in  every  case  an  equation  of  the  form 
By  Y  =  y  +  A  must  be  solved  f or  y .  This  equation  is  encountered  also 
in  the  first  doorway  if  Method  F  is  used.  In  Table  E-2  the  meanings  of 
y,  A,  and  B  for  each  case  are  listed.  Since  conditions  in  the  first  door 
way  under  Method  F  are  independent  of  conditions  inside,  the  definitions 
for  the  first  duct  appear  separately. 


In  Figure  E-11  results  of  a  calculation  carried  out  by  Method  F  for 
two  rooms  are  compared  with  measurement®^  of  average  pressure  in  the  firs 
room.  The  two  rooms  consisted  of  two  small  models  such  as  those  illus¬ 
trated  in  Figures  E-5  and  E-6  placed  back  to  back  with  a  connecting  door 
exactly  like  the  outside  door.  The  experiment  was  performed  in  a  shock 
tube  in  which  the  wave  struck  the  first  doorway  head-on.  Also  in  Figure 
E-11  appears  the  result  of  a  calculation  treating  the  whole  volume  of  the 
two  model  rooms  as  if  it  were  in  a  single  room.  In  Figure  E-12  measured 
and  calculated  results  for  the  second  room  in  the  model  appear.  At  least 
in  this  one  example,  all  three  calculated  pressure  rises,  i.e.,  in  the 
first  room,  in  the-  second  room,  and  in  the  whole  volume  of  both  rooms 
treated  as  a  single  room,  are  quite  similar  and  there  appears  to  be  no 
advantage  in  using  the  complicated  procedures  for  computing  the  fill  of 
two  connected  rooms . 


The  calculated  histories  shown  in  Figures  ii  and  12  have  not  been 

carried  beyond  the  time  of  equilibrium  between  inside  and  outside  pres- 
sures . 


Figure  E-13  shows  pressure  history  calculated  by  Method  D  ignoring 
the  wall  between  rooms.  For  this  calculation  the  discharge  coefficient 
has  been  set  equal  to  0.7  on  inflow  and  1.0  for  outflow.* 


CAVFIL,  a  FORTRAN  program  written  at  IIT  Research  Institute  was 
used  to  make  the  computation  ^ 
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Table  E-2  (concluded) 
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FIGURE  E-12  COMPARISON  OF  OBSERVATIONS  IN  THE  SECOND  ROOM  OF  A  TWO' 
ROOM  MODEL  WITH  CALCULATIONS  BY  METHOD  F 
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FIGURE  E-13  COMPARISON  OF  OBSERVATIONS  IN  TWO-ROOM  MODEL  WITH 
CALCULATIONS  BY  METHOD  D 
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Openings  into  Different  Pre 


ssure  Fields 


The  simple  methods  considered  unriev  i-ur.  i,  ,  • 
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A.  Computer  Program 
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FIGURE  E-14  SKETCH  ILLUSTRATING  NUMERICAL  EXAMPLE 

head-on  incidence  upon  the  opening  to  the  larger  room,  and  positive  phase 
duration  of  1  sec.  Standard  ambient  conditions  will  be  assumed,  i.e., 

Po  =  14.7  psia  and  =  0.0761  Ib/ft^.  In  the  front  wall  the  opening 
area  is  25  ft  and  there  is  a  like  opening  in  the  rear  wall.  Because  the 
door  between  the  rooms  has  an  area  comparable  to  the  area  of  each  of  the 
outside  doors,  the  presence  of  the  partition  will  be  ignored  and  the 
volume  to  be  filled  taken  as 

3 

=  50  X  20  X  8  =  8000  ft 

As  a  first  step  the  pressure  histories  outside  the  two  doors  will 
be  calculated,  according  to  the  procedures  recommended  by  Ref.  1.  Since 
the  first  window  is  struck  head— on,  there  will  be  no  time  delay  there 
and  the  peak  pressure  will  be  the  reflected  pressure  P„  which  is  calculated 
as : 


K-67 


Here 


IV  +  '1  I' 
o _ 

P_  =  2  *  7  ^  +  p 

R  o  so 


p  =10  psi  and  P 
so 


=  14.7  psi,  hence 


p  =25.3  psi 
R 


.e  will  be  felt  in  declining  strength  for  Iho  duru- 
This  reflected  pressuie  nv  b  ^  estimated  as 

tion  of  the  clearing  time,  t^,  which 


3s 

"^c  =  T 

-  r>f  the  wall  undergoing  pressure 
where  s  is  a  dimension  of  the  wax 

sound  speed.  In  this  case 


clearing  and  c  is 


-  [v  P  /p 


144/0.076096]^''  =1118 


and  s 


so  that 


t  =53.7  ms 
c 


47-.  uressure  during  the  interval 

decline  of  reflected  pressure- 
As  an  approximation  the  decli 

0  ^  t  ^  t^ 

.  t  -  t  the  pressure  Pe  on  the  first  wall 

is  treated  as.  linear,  that  overpressure  P_  and  the  drag  P^^ 

is  simply  the  sum  of  the  ree  ^  the  time  ^ 

arising  from  the  exponential!/ from  P^^  =  10  Psi 

free  field  pressure  has  fallen 


P  =  P 
sc  so 


t,  \ 
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where  is  the  free  field  duration  of  positive  overpressure;  in  this 


case  t^  =  1  sec  so  that 


P  =  8.97  psi 
sc 


Peak  drag  pressure  P  is  computed  from  the  formula 

do 


u 


do  2 


7  P  4-  P 
o  so 


5 

2 


100 


7  X  14.7  +  10 


=  2,21  psi 


so  that  drag  pressure  at  t  =  t  is:^ 

c 


p  .=  p  I  1  .  ^  r  e(-2 

dc  do  I  t 

o 


Numerically  this  is 


P  =  1,78  psi 
dc 


Therefore  the  pressure  outside  the  first  opening  at  t  =  t  is; 

c 


P  =  p  +  p  =  8.97  +  1.78  =  10.75 
c  sc  dc 


and,  assuming  a  linear  fall  from  Pg^  to  P^ ,  for  0  £  t  ^  t^ ,  pressure  P^^ 
outside  first  wall  as  a  function  of  time  becomes; 


(t  -  t\ 

P  =  P  +  ^ - -  (p  -  P 

1  c  t  \  R  c} 


(68) 


For  the  remainder  of  the  positive  phase  duration  outside  pressure  at 
the  first  opening  is  simply  the  sum  of  the  decayed  side-on  pressure 


P  =  P 
s  so 


,  t  \  -t/t 
1  -  —  \  e  o 

^  I 

<  o  / 


(69) 


and  decayed  dynamic  pressure 
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so  that  for  t  >  t  >  t 
o  c 


p  =  p  +  p 
1  s  d 


(70) 


(71) 


For  the  first  wall  a  drag  coefficient  equal  to  +  1.0  has  been  tacitly 
assumed  above  but  for  the  second  opening  this  coefficient  will  be  different 
from  1,  and  according  to  Ref.  1,  a  value  of  -0.4  may  bo  assumed.  At  the 
rear  opening  there  is  no  reflected  pressure  and  the  delay  is  equal  to  the 
time  taken  by  the  front  to  traverse  the  building  (assuming  that  the  open¬ 
ing  is  already  open  upon  blast  arrival  or  that  it  is  immediately  forced 
open  by  the  blast) .  Blast  front  speed  can  be  found  from  the  formula 


1  + 


Y  +  1 
2 


wh0r©  c  =  sound  sp0©d  in  Eitibi0n't  air ,  i  #0  •  f 
o 


c 

o 


Y  P  32  X  i44/p^ 


1/2 


=  1116  ft/s0C 


H0nc0 


U  =  1116  X 


2.4  10 

2.8  ^  14,7 


1/2 

=  1148  ft/s0C 


so  that  delay  at  the  rear  entrance  is  =  43.5  ms 

Beginning  at  t  =  43.5  ms  the  room  starts  to  fill  through  the  rear  opening. 
(Outflow  through  openings  other  than  the  first  can  ordinarily  be  neglected 
during  the  delay  period:  either  the  other  openings  are  closed  to  the  blast 
for  a  certain  period  or,  if  not,  the  blast  travel  time  to  them  is  much 
shorter  than  the  time  required  to  start  outflow  through  them.)  Filling 
through  the  rear  opening  takes  place  however  from  a  reservoir  at  lower 
pressure  than  that  outside  the  front  opening,  i.e.,  outside  pressure 
at  the  rear  is: 

p=P— 0.4P  (72) 

Ir  s  d 
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The  decline  in  Pg^  and  which  occurs  while  the  blast  front  travels 
from  front  to  rear  opening  is  negligibly  small  and  may  be  safely  neglected 
for  buildings  of  ordinary  size. 

For  this  sample  case,  the  quantities  P^  (pressure  outside  the  front 
opening)  and  Pj^^.  (pressure  outside  the  rear  opening)  have  been  calculated 
as  functions  of  time  and  plotted  in  Figure  E-15.  The  figuro  shows  the 
discontinuity  in  the  derivative  of  Pj  with  respect  to  time  at  t/;.<; 

(^c>  when  the  reflected  pressure  is  assumed  to  disappear  and  the  out¬ 
side  pressure  takes  on  its  quasi-steady  value. 

Also  plotted  in  the  figure  are  inside  pressure  histories  P3  calculated 
by  two  methods:  the  greatly  simplified  procedure  given  in  section  II-A 
of  this  report  and  the  step-by-step  method  F  explained  in  Section  II-B. 
According  to  the  first  of  those  two  pi’ocodures  the  estimated  history  is 
given  by  the  lino  segments  ODFG  obtained  as  follows.  \Vhen  the  blast  ar- 
v'ivos  at  the  front  opening,  filling  immediately  begins  along  line  OA, 
vhorc  point  A  is  the  intersection  of  the  outside  pressure  history  and  the 
abscissa 


V 

3  8000 


2A^  2  X  25 


160  ms 


Line  BC  is  a  similar  line  representing  filling  through  the  rear  opening, 
beginning  after  the  delay  time  of  43.5  ms.  Ordinates  under  the  line  BC 
have  been  added  to  ordinates  of  OA  to  produce  the  line  DE.  Since  areas 
of  both  openings  are  equal,  the  point  F  is  placed  halfway  between  current 
outside  pressures  outside  the  front  and  rear  openings  and  the  decline  of 
P3  represented  schematically  by  the  line  FG. 


The  step-by-step  calculation  results  in  the  curve  labelled  "Method  F" 
in  Figure  15.  Because  of  the  high  reflected  pressure  during  the  interval 
0  t  ^  tc  (which  does  not  influence  the  results  of  the  simplified  method 
in  this  example)  the  more  careful  calculation  shows  a  faster  build-up  of 
room  pressure  than  the  line  ODF.  To  demonstrate  the  method  the  first  step 
of  the  stepwise  solution  will  be  calculated  below. 

Since  the  least  sound  transit  time  across  the  room  is  approximately 
20  ms  we  will  choose  a  value  of  At  =  5  ms.  At  t  =  0  there  is  only  one 
opening,  that  in  the  front  wall.  Outside  pressure  there  at  that  time  is 
Pr  =25.3  psig.  Inside  pressure  P3  =  P  =  14.7  psia  and  density  is 
p'  =  p  =  0.076096  lb/ft3.  ° 

*3  O 
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1. 


=  25,3  +  14.7  =  40.0  psia 


P  l/v  r  1 

o  1  40.0  1/1.4  ^  : 

2.  P  =  —  p  =  (0.076096)  =  0.156  Ib/ft 


3.  P^  =  0.1912  P^  =  (0.1912)  (40.0)  =  7.65  psia 


P2  =  F 


1^2  1/Y 


=  (0.307)  •  (0.156)  =  0.04785  Ib/ft" 


2  2 

5.  u  =  — 
2  V  - 


P  P. 


-  ijl 

o 

b 

7.65 

“  0.4  1 

0.156 

"  0.04785 

32  X  144  =  3.11  X  10 


6  ft 


1  2  rr  .r.  Ib/ft  . 

6.  q^,  =  “  u  p„  =  7.45  x  10  - x-  (dynamic  pressure  at  front  opening) 

21  2  2  2  sec^ 

7.  Since  P,  >  P^:  u  >  0,  i.e.,  flow  is  inward.  Were  P^  >  P  ,  u  would 

1  3  2  3  1  2 

be  negative 

8.  Am  =  u  p  A  At  =  0.0105  lb 

0 1  ^  ^  td 


1  *^31  (40.0)  (0.0105) 

9.  Aw^,  =  7 - -T -  =  - - -  •  144  =  1360  ft  lb 

31  /  V  -  1\  0.  0.4  ... 

V")  ^ 


10,  Since  the  rear  opening  is  closed  at  this  time: 


Am  =  Aw  =0 
32  32 


W'ere  the  second  opening  available,  steps  1-9  would  be  repeated 

using  initial  outside  pressure  at  rear  opening  to  calculate  Am 

and  Aw  ^ 

32 

11.  Am  =  Am  +  Am  =  0,0105  lb 

*3  31  32 

12.  Aw  =  Aw  +  Aw  =  1360  ft  lb 

3  31  . 32 
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13. 


/  (Y  ” 

=  p'  +  ' _ 


Aw- 


3  3 

pressure) 


14.7  + 


(0.4) (1360) 
8000 


14.768  psia  (room 


in  room) 


0.076096 


9  =  0.0760973  Ib/ft^ 

8000 


(air  density 


15. 


Set  equal  to  P^ 

and  p'  equal  to  p  and  return  to  step  1  with  value  of 
P  at  t  =  At  =  5  ms  • 


V.  Edge  Diffraction  of  an  Acoustic  Wave 

A  weak  planar  shock  striking  a  semi-infinite  wall  head-on  can  be 
treated  approximately  as  a  self-similar  acoustic  wave  in  the  manner  demon, 
strated  by  Ludloff.^®  Such  treatment  is  two-dimensional  and  neglects  the 
presence  of  floor  and  ceiling. 

In  the  acoustic  approximation  all  disturbances  or  effects  are  propa 
gated  with  sound  speed.  Thus  after  the  incoming  wave  front  strikes  the 
semi-infinite  wall,  the  influence  of  the  edge  will  be  felt  only  within  a 
cylinder  whose  axis  is  the  edge  and  whose  radius  is 

c  t 


where  c  is  sound  speed  and  t  is  elapsed  time  after  initial  impact. 

If  a  cartesian  coordinate  system  is  placed  so  that  the  edge  becomes  the 
z-axis  and  the  negative  y-axis  lies  in  the  wall  the  location  of  the  circle 
of  disturbance  due  to  the  wall  is 


2  2  2  2 
X  +  y  =  c  t 


and  the  equation  satisfied  by  the  overpressure  p  is 


2  ,  2 

5  P  1  5  P 


2 

5y 


2  ^2 

c 


If  a  change  of  variables  is  made: 


(73) 
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X 


(7^) 


ct 


and  tan  G 


z  =  H 

X  c 


then  in  the  new  coordinates  (T|,  a,  6)  the  di stiirl^anco  in  confined  to  n 
circle  of  unit  radius.  If,  further,  radii  are  changed  in  scale  according 
to  the  formula: 


P  = 


1  + 


<1  - 


(75) 


where  p  is  the  new  radius  and  r  =  T]  +  a  ,  then  the  equation  satisfied 
by  p  in  the  cylindrical  polar  coordinates  (p ,  9)  can  be  written 


(76) 


This  is  Laplace *s  equation  and  is  satisfied  by  the  imaginary  component  of 
any  analytic  function  of  the  complex  variable. 


In  the  case  of  edge  diffraction  of  a  weak  shock  the  present  method  yields 
a  solution  for  0  r  ^  1  (or  0  ^  p  ^  1)  ,  The  angle  9  lies  in  the  range 

-  to  as  illustrated  in  Figure  E-16.  The  origin  of  coordinates  is 

A  A 

labelled  0. 


The  boundary  conditions  are  determined  by  physical  considerations. 
Because  the  acoustic  Eq.  (73)  or  (76)  is  linear  in  p,  the  incident  over¬ 
pressure  may  be  taken  as  unity  and  the  pressure  reflected  from  the  wall, 
as  2.  The  overpressure  in  the  undisturbed  air  is  of  course  zero.  These 
conditions  imply  the  following  pressure  values  on  the  circumference  of  the 
circle  of  disturbance 


J  ^  0  ^  0 
2 


0  ^  9  TT 

^  A  Srr 
TT  ^  0  ^  — 


p  =  0 

p  =  1 

p  =  2 


(77) 


The  areas  of  uniform  pressure  are  marked  in  Figure  E-16, 
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{-PLANE 


INCIDENT  WAVEFRONT 


P  =  l 


P  =  0 


REFLECTED 

WAVEFRONT 


(UNDISTURBED  AIR) 


FIGURE  E-16  BOUNDARY  CONDITIONS  ON  CIRCLE  OF  INFLUENCE 


Potential  theory  guarantees  that  there  is  only  one  pressure  distribution 
within  the  unit  circle  which  satisfies  Eq.  (76)  and  meets  the  foregoing 
boundary  conditions.  Solution  of  Eq ,  (76)  will  be  guaranteed  by  the  fact 
that  the  distribution  function  is  analytic. 


A  further  transformation  of  variable  will  make  the  exposition  of  the 
correct  distribution  function  clear.  Let 


r,  1/2 

where  R  =  p 


and 


$ 


(78) 


In  the  w-plane  the  unit  circle  of  disturbance  becomes  a  semicircle  of 
radius  1;  the  back  side  of  the  wall  lies  along  the  line  §  =  0  while  the 
front  side  falls  along  §  =  rr.  The  line  0  =  0  in  the  £-plane  rotates  to 
$  =  ^  and  the  line  0  =  n  is  found  along  §  =  Thus  in  the  w-plane  the 

boundary  conditions  on  the  semicircle  R  =  1  become 


E-76 


p  =  0 


TT 

0  $  <:  - 
4 


p  =  1 


^  ^  S  ^  3tt 
4  4 


p  =  2 


^  ^  n 


Any  analytic  function  of  w  will  also  be  an  analytic  functi 


;  t.ion  ot 


The  function  w  -  e  is  represented  by  the  line  BA  in  Figure  e-17 
and  its  argument  by  the  angle  <y.  imen  A  falls  on  the  unit  circle,  a  in- 

TT  3  TT 

creases  discontinuously  from  -  ~  ^  moves  counter  clockwise  through 

the  point  B.  The  function  w  -  e  ^4  is  continuous  within  and  on  the  upper 
semicircle  DBE.  Furthermore  the  included  angle  y  is 

Y  =  O'  +  P 


where  Arg  w  -  e 


O'  and  Arg  [w  -  e  J  =  P.  But  along  the  arc 


Y  =  7*  Hence  the  function 


^  ^  /  A  w  -  e 

f ,  =  -  \  Arg  - 

1  TT  TT 

1— 

,  4 

'  •—  w  -  e 


is  zero  along  EB  and  1  along  BD. 

By  a  similar  argument  based  on  Figure  E-18  it  is  clear  that  the 
function 


^  /  A  w  -  e  tt 

f  =  —  <  Arg  - - 

2  TT  .3tt  2 

^4 

'  w  -  e 


is  zero  along  the  arc  Eb'  and  1  along  b'd. 
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Thus  the  sum 


meets  the  required  boundary  conditions  along  the  semicircle  EBD  in  the 
w-plane .  But 


(.5-n  .7tt 

.  tt  .3ti 

1—  1 — 

4  4 

•-  w  -  e  w  -  e 


is  analytic  in  w  and  z  and  since  in  Eq ,  (80)  is  evidently  the 

imaginary  part  of  F(w) ,  Eq .  (80)  provides  the  sought-for  expression  giving 
the  px’ossure  distribution  within  and  on  the  circle  of  disturbance  p  =  R  =  1, 


Eq.  (80)  can  be  written  as 


TTp  =  tan 


“  *  ^/5  \  -1  r  ^  *  :/l 

- ‘Y  ■  I  +  tan  I  - ^ 

R  cos  $  -  /r  /  \  R  cos  §  - 

^/2  /  \  V2 


“  tan 


R  sin  $  - 

J2 


R  cos  $  - 


“  tan 


-1  V  ®  “  /2 


R  cos  §  + 


where 


R=  p 


1  .  (1  - 


TT  0 

^  =  T  +  t: 

4  ? 

The  denominators  of  the  four  arctangents  in  (81)  have  I'eal  roots  in  the 
range  0  ^  ^  ^  n  when  p  s  1/2;  and  care  must  be  taken  to  evaluate  the  in¬ 
verse  functions  in  such  a  way  that  p  is  continuous  within  the  unit  circle, 


There  are  no  zeros  in  the  denominators  when  o  <  “  • 

^  2 
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then 


and 


If  R  cos  ^ 


B  cos  ^ 


3rr  .  ^ 

—  ^  ^  tt 

4  1 


rr 


2  4 


0  and 


0 


If  the  FORTRAN  algorithm  is  used  to  evaluate  the  arctangent  terms  in  (81) 
then  for 


the  second  term, 


$  >  §  and  R  > 
2 


^/5 


tan 


.  R  sin  §  +  ,- 

-1  /  _ J2 

1 

R  cos  §  - 

^2 


may 


be  increased  by  rr,  and  the  third  term 


tan 


-1  I  "  -  ./5 


R  cos  ^  ^  r 

V2 


must  be  decreased  by  tt,  or  the  sum  in  (81)  increased  by  2tt.  \V\\cn 
§  >  and  R  ^  another  2tt  to  keep  the  expression  f< 
continuous  and  insure  the  existence  of  its  derivatives. 


In  order  to  assure  continuity  with  respect  to  radius  the  foregoing 
choices  limit  the  arctangents  to  certain  quadrants  where  radius  R  is  such 

that  no  zeros  in  the  denominators  exist,  i.e.,  when  ^  Thus,  since 

we  have  chosen  to  add  it  to 


-1 

tan 


R  sin  5  +  ^2 


R  cos  §  - 


^/2 


T 


foi’  $  >  and  R  >  “,  we  must,  when  using  the  FORTRAN  algorithm,  add  tt 
to  the  same  inverse  tangent  for  all  §  when  R  <  i  in  order  to  place  the 
^oiuputccJ  in  'tl'10  sGcond  QUcidi^ciri't ,  Sirnilcii*ly  ^  by  subir’ctcii  S 


when  §  >  §  and  R  >  -  we  are  placing  the  arctangent  in  the  range  -  -  to 

^  ^  o 

1 

-  rr;  hence  the  same  function,  when  R  <  “>  computed  in  the  same 

range  by  subtracting  n  from  the  angle  as  computed  by  FORTAN.  When  R  <  - 
the  denominator  ^ 


R  cos  $  +  r 

V2 


is  always  positive  and  when  R  <  ^  the  FORTRAN  choices  of  quadrant  for  the  two 

arctangents  containing  it  are  consistent  with  those  for  R  >  — . 

2 

Pressure  contours,  located  as  outlined  above,  are  shown  in  Figure  E-19. 
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NOTATION* 


A  area  of  opening  or  openings  through  which  room  fills 

A  constant  equal  to  various  functions  of  y,  appearing  in 

Equation  (5) 

B  constant  equal  to  various  functions  of  y  and/or  of  interior 
and  exterior  conditions,  appearing  in  Equation  (5) 

c  speed  of  sound  in  air 

Cy  specific  heat  at  constant  volume  for  air 

F  numerical  factor  for  consistency  of  units 

I  integral 

K  discharge  coefficient 

distance  along  jet  axis  measured  from  opening 
m  mass  of  air 

Am  increment  of  air  mass  corresponding  to  time  increment  At 

P  air  pressure 

air  pressure  during  immediately  preceding  time  step 
^lo  P®ak  overpressure  outside  front  or  first  wall 
total  overpressure  outside  rear  wall 
Pc  sum  of  dynamic  and  side-on  air  pressure  at  time  t=  tc 
^d  free-field  dynamic  pressure 

*  Numerical  subscripts  may  be  used  with  variables  to  refer  to  a 
particular  space  or  volume;  see  Preface, 


^dc 

dynam: 

jessure  i' 

a  St  wave  at  time 

^do 

peak  d 

lie  presr 

in  blast  wave 

Po 

ambien 

’  pressr 

Pr 

ref  led 

res sure 

1  wall  struck  last  wave 

Ps 

f  ree- 

side-oi 

pressure 

^sc 

f  re^ 

side- 

.>verpressur  ..  t=  tc 

^so 

free-: 

peak  B 

verpressurc 

p 

acoust: 

ver) prr 

q 

dynamic 

.sure  c 

q 

(spatii 

/erage 

c  pressur  ,  . 1  at  one  cross  section 

R 

gas  co- 

c  for  r 

:e  Ref.  5.) 

R 

radius 

)t  of  £ 

ng  air 

r 

radial 

linate, 

.  r  =  T12  +  c; 

r 

radial 

iinate 

od  from  axis 

S 

specif J 

i:ropy  c 

AS 

increm< 

f  entrr 

s 

wall  d 

■>ion  use 

'ompute  clearin. 

T 

absolu 

emperatui 

AT 

fillii 

me  or  int 

between  firs  al  of  blast  and 

achic”' 

:t  of  prcK' 

qui librium 

t 

time 

ired  from  : 

arrival  of  it  structure 

clearing 

time  of  structure  in  blast  wave 

duration 

of  positive 

side-on  overpressure 
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At  increment  of  time 
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time  in- pendent  coordinate,  .c.,  'fl  -  y/ct 

1  ^-'iirface  and  the  positive  direction 
angle  Ik  tween  inward  normal  .  hurtace  ana  f 

of  the  x-axis 


angular  coord 


inate,  i.e.  ,  tar-  -  y/x 


density  of  air 


reduced  radial  coordinate 


air  density  during  immediately  preceding  time  step 


ambient  air  density 


time  Independent  coordinate,  i.s.,  <1  -  *'<=1 

time  required  to  tranemlt  a  pound  alsnal  over  the  longept  room 


dimension 


angle  in  complex  w-plane 


! 

f 
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